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HY DRAULICS DIVISION 


of the American. Society of Civil ngincers 


ENGINEERING APPRAISAL OF HYDROLOGIC DATA 


Progress Report of the Task Group 
Hydrologic Data of the Committee on 
Hydrology of the Hydraulics Division 


is is given of ‘studies 1 made 1936 on the adequacy 


= 


and engineering groups reveal a Se need for more and better il 


‘hydrologic d data. The resulting reports by governmental, non-governmental 


_ processing and need to accommodate new types of data. ae 
7 ‘polut to strengthening data programs at both state and federal levels. The 


responsibilities of civil engineers are outlined. 


age: stony, 


=p Historical Bac 
‘Task Group was established by the Committee on in late 
1955 for the purpose of obtaining an objective appraisal, from an engineerin 
_ viewpoint, , of the adequacy of hydrologic data and the current programs of | 
collection with respect to present and anticipated needs. The Task Group 
_ held one meeting, in October 1957, at which preliminary drafts of of the e report 
reviewed and an outline of the final report was approved. 
_ Several studies have been made in recent years in which the adequacy of 
the nation’s hydrologic information has been examined. With two possible _ 
_ exceptions (1949 report of Federal Interagency River Basin Committee and | 
_ 1951 report of Engineers Joint Council), most of these studies treated hydro- ; 
4 logic data only incidentally, or as part of other and perhaps la larger problems. = 


Policy or the “ certain federal functions. 


a of the copyrighted Journal of the Hydraulics Division, Proceedings of 
American of Civil Engineers, Vol. 85, No. 7, July, 1959, 
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= useful and valid. Others pointed out certain deficiencies, many of which 2 
been overcome in recent years. All problems o of mee adequate — 
a have mot, of been solved. 


This report should be useful to the Society member and to. 
neering profession in general because: 


a It is a current engineering appraisal of hydrologic | data in in the United 


x 2. It contains a key to the sources of available hydrologic data. "a eee 
It reveals recent improvements in federal data collection programs and 
further improvements which may ‘result from current reappraisals. 
_ 4, It points out opportunities for making better use of available data and | 


adding to the volume of public records, | 


5. It reiterates the civil engineer’s responsibility for 


rs data are used in nearly all walks of life and for many varied me 


of hydrologic data to him is indicated a definition of civil. engi- 
' 6 neering taken from the California statute regarding registration of engineers; a 


Civil engineering embraces the investigation of the laws, phenomena an and 
forces of nature in connection with fixed works for irrigation, drainage, 
_ Water, power, water supply, flood control, inland atatittas harbors, 
‘Iti is obvious that the requirements for hydrologic and the 
ity for its collection a and analysis area problem of sufficient it importance t to ~ 


merit close attention by the civil engineering profession. 
Status of Hydrolog ic Da ta a 


ee From the point of view of the practicing civil engineer, hydrology is tbe 
science of the occurrence, movement, and control of water. _ The engineer: is 
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.* Hydrologic ron consist, therefore, of the knowledge gained through « obser 

pt vation and study of the | phenomena in the hydrologic cycle. These become the | 
tools of the hydrologist and the engineer. ed into the 


categories: 
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are observed data which have been converted or otherwise put into usable or 


ae Analyzed data include basic data which have been subjected to elemen- 


tary degrees of engineering treatment within the limits of the curve. obser- 


vations. Examples are map and the curve 


: sult of interpretation and extrapolation beyond the limits of observation, 
through application of sound engineering procedures. An example is the esti- — 
ee yield of ar an aquifer. Such data may have ) many applications | to engineer 


in blem 

Basic data, however, usually are much more than mere statistics. Some _ 

+ _ experience a and skill goes 3 into” even the most rudimental observations. . Mere oa 

statistical records often are valueless and possibly misleading unless knowl- 

; edge of outside influences has been applied. = ‘Thus, most observational data = 


engineering g problems. The ‘following are considered to be the principal types 
stream flow sediment load of streams 
river and lake levels eer 
Closely related to hydrologic data are other observations and technical “a 
- knowledge in the related fields of meteorology, agronomy, geology, and — ; 


~ 


sciences. factors are useful to engineers i in 


Present Status of Hydrologic Data 


scope of water problems and with subsequent modifications in programs 


a water investigation. Thousands of items of hydrologic data are being ob- 
a served, collected, and developed each day throughout the entire field * aaa 
hydrology. iItis estimated, for example, that more than three million Roms 
v7 of basic > hydrologic data 2 are collected each year i in n California alone. ‘The © 


sound engineering principles to the use of data. 
_ n fulfilling his duties as a director of hydrologic data programs, the en- 
_ gineer may be required to broaden his vision to encompass both a rapidly 
_ growing science and a new concept of its importance. - While the collection ate 
of data at present is generally for a specific purpose, 1 we cannot today ac eo 
curately predict the use or purpose to which today’s data may be Coe in % 
The present status of hydrologic data isa remarkable 
_ sidering that the United States was already in its third century of history a 
_ when systematic public collection of river stage and rainfall records beg 
_ Relatively speaking, however, most of the records are short- term. —" | 


Basic data, generally speaking, are observed values or measurements. q yal 
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of gages for measuring stream and precipitation now constitute 
- both a systematic and a long-range approach to the collection of hydrologic .: - 
data. With the growing awareness that data needs are constantly changing, ie. 
= networks themselves should receive periodic reappraisals. 

_ Research also contributes to the knowledge of the physical processes in- | J 
“volved in the water cycle. . This knowledge helps to determine what features” an 
- should be measured in the basic data programs. The entire field, therefore, 
by as never before in the short of hydrology. ( (1) 


‘Present Status of Collection Programs 2 


Hydrologic data are ¢ collected by Federal, state, and 


agencies of government, and by private firms and individuals. The 
Federal Inter-Agency Committee on Water Resources has compiled a 
d 


= 


plete record of the federal programs in a report which has had only limite 


distribution. The summary is based largely on that report. (2) 


’ a Within the Department of Agriculture, basic hydrologic data are collected 
i in the Agricultural Research Service, the Soil Conservation Service, and the — 


| - Forest Service. The data are needed to em, on the authorized functions of 


The Agricultural Research Service collects basic data on precipitation, ag 
runoff, infiltration, ground water, erosion, sediment yield and reservoir ‘silt-— 
ing, soil r moisture, evaporation, and consumptive u use, and makes interpreta aa 
tions of these data for solution of agricultural problems. wi 

_ - The Soil Conservation Service collects basic data on irrigated lands and, in 

portions of the West, sponsors and coordinates snow surveys. also collects 

be _ some data on consumptive use, floods, and infiltration for use ‘in project de - 

sign. _ The SCS sometimes allocates funds to the Geological Survey and the od 
‘The Forest Service collects basic data on precipitation, interception, infil- ame 
~ , runoff, humidity, and consumptive | use at experimental forests : and 
= ‘It collects data on snow, floods, and reservoir sedimentation ” no 
forested areas for use in good forestry practices. 
@ _ The Corps of Engineers, U. S. Army collects hydrologic data for igs 
- and operating its flood control projects and navigation facilities. In addition, — 
o it allocates over a million dollars annually to the Weather Bureau, the Geologi- 
cal Survey, and other agencies for the collection of basic data on precipitation, 
—— sediment load, and ground water. The Corps collect basic data .- 


7 on precipitation and streamflow at sites where the data are needed in current — : 
operation; data on wind, tides, and evaporation for use in design; stage data _ 
_ at locks and dams; and ‘sediment - load records on streams flowing into iil 


q 
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2 reservoirs. As part of its work for the Mississippi River Commission, the 


data on hydrologic factors the level of the Great Lakes. 
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artment of Commerce | 


Within the Department of , the Coast and Geodetic Survey 


the Weather Bureau are authorized to collect and publish hydrologic data for 
others as well as primarily for theirownuse. 
| ‘The Coast and Geodetic Survey collects basic data on tides to predict their — 
‘time and height in the tidal reaches of major rivers. = 
&s The Weather Bureau is one of the major collectors of k basic hydrologic data. 


‘they a are all used widely in hydrologic analysis. The 


soil-moisture conditions, particularly as these affect 


ture. it has responsibility for forecasting | flood stages and for ‘making water-— 


supply forecasts in the West and in the upper Missouri River basin. 


pollution. It reports on waste disposal, municipal use of water, biological, d 
7 “bacteriological, and chemical quality of surface and groune water when as- 


Within the ‘Department of basic data are collected b 
; the Fish and Wildlife Service and the Bureau of Reclamation for ve mae aq 
“operating their projects, and by ‘the Geological Survey. 
_ The Fish and Wildlife Service measures the temperature of surface water 
for use in connection with fish-hatchery 
= The Bureau of Reclamation allots ‘money to basic data agencies such as the 
- Weather Bureau and Geological Survey for collecting data on precipitation, ae 
evaporation, streamflow, ‘sediment load, and ground water for existing or pro- 
- posed reclamation projects. The Bureau of Reclamation collects data on con- 
sumptive use, storm precipitation, and historic floods for designing and oa 
operating irrigation and power projects in the AT western States. ‘It also col- 


lects discharge records on some of its canals; precipitation records at its 


4 projects; and, in cooperation with other agencies, basic data on density cur- — 
reservoir evaporation and sedimentation, and snow melt. 
‘The Geological Survey, another basic data agency, collects, , analyzes, and i. 
 pebten many types of hydrologic data for others to use in developing water _ 
‘supplies, abating pollution, and protecting against floods. It is the principal _ a 
af i source of basic data such as stream flow, sediment load, chemical quality, _ 
and ground water. It reports on floods and droughts as they affect ne 
and ground water. It collects data on water temperature, consumptive use, 
and reservoir sedimentation, and collaborates with other agencies in collect-_ 
ing data on density currents, evaporation, snow surveys, and tides. 
aa es Tennessee Valley Authority ¢ collects hydrologic data in connection with ; 
the the operation of reservoirs and area | development. It It collects and publishes 
. a data on: 1 reservoir stages, evaporation, density currents and sedimentation; a 


infiltration, interception, soil temperature, soil moisture a and ground water at 
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voirs and ine. “Much of the stream gaging in the TVA area is financed a : 
transfer of funds to the Geological Survey. 
_ The International Boundary and Water Commission collects hydrologic data 
connection with carrying out the provisions of the Mexican treaty for ap- 
propriation 1 of the water of the Rio Grande. It measures streamflow, | floods, — 
sediment load, chemical quality, and reservoir storage on the Rio Grande and 
tributary streams between Elephant Butte Dam, New Mexico, and the Gulf ae 
_ Mexico. In addition, it collects some data on precipitation, air temperat ure, a 


a humidity, and wind in that area. In me Lower Colorado River basin it collec 


State geological surveys, weter engineers and 


a - other agencies of state government either cooperate with federal agencies in 


‘ti the cooperating parties share both in the financing and in the engineer- 
_- work, while in others the cooperation is entirely financial. Most of the 
; states having cooperative programs with federal agencies rely on the federal 
ape publications to release the statistical data and confine their publications to 
summaries or interpretive reports. Some state agencies collect data which > 
is of state interest only, such as might be obtained in a census. States often | ; 
: _ obtain data through research and experimentation at state universities. a : 
Additional data of state interest are gathered by state organizations having © 
legal or administrative responsibilities in conjunction with water problems. 
of | meridian that adminis- 


exist for a number of major in the United 
States. _ Generally the of data is done through the compo- 


: Regional and and Local 


special districts, etc. Some of these agencies basic hydrologic 

required for the ‘fundamental purpose of the organization. For example, sani- 
_ tary districts collect data on stream height and flow, waste water a 
- water quality, and weather. Such local data may be exceedingly valuable to _ | 


. carrying out regional hydrologic studies, but it usually is on 1 file only red 


Water supply authorities, districts, and utilities serving municipalities as- 
semble data on ground water and surface water sources being used. They | 
a seldom assemble data on sources not being used, and in many instances data 
are lacking even on the sources bein used. 
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Private Firms and Individuals 


ee may collect and maintain records of hydrologic data for their — 
+ use. In nearly all cases the collection is related to the specific business 


activity and is therefore specialized both as to type and location. 


_ Engineers and consultants in private practice are not ordinarily silane 
: of basic data for public use. Some engineers, however, make their records “2 
available to collecting agencies, but few are published. Private | records and 
data from public sources are used by engineers and consultants | to develop 
‘new or related data. ‘Such data usually become the property of the client and 
are never published. may appear in technical articles 


Policies Data Collection 
* eau Flowing water differs from land | and o other fixed property in the important 
4 ‘respect that it cannot be retained for use by o one } person but passes on to ‘other | ; 
users, thereby creating complex problems in water rights and damages. vas This a 


The Geological Survey the Weather: Bureau two agencies 

_ primary federal al responsibility t to collect and publish hydrologic data. The - 
_ policies which govern their operations have been described in statements = 


submitted by their chief executives to the Task Group. state- 


ments are included in the © Appendix. pf 


Authority 


OLS 


Geological Survey was in the of the Interior in 


1879, but its responsibilities Soe. wakes resources investigations were not de- Ag 


until several years later. The first appropriation was made in 1888, — 
_ and successive appropriation acts have continued to define its authority i, 
- “gauge the streams and determine the water supply of the United States, in-— Bi 
4 cluding the investigation of under-ground currents and artesian wells in arid 
and semi-arid sections. ” Its cooperative investigations with the states have 
been a major part of the work for more than 50 years, and wupeaia on this a 
Objectives of the Geological Survey’s programs are the determination and 
description of the quantity and quality of the surface and underground waters Mi 
int of the country — an appraisal of the total water resources. Over the years Load 
the emphasis has shifted from mere collection of basic data to the nos somal % 
and interpretive phases of water resources investigations, ‘Research and ex- 
_ perimentation are becoming important adjuncts of the Survey’ s investigative 
responsibilities. ultimate objective is an improved knowledge of natural 
phenomena as the basis for planning and design of water resources projects. © 
_ The Weather Bureau was organized under the Department of Agriculture rn 


. E in 1890 and was transferred to the Department of Commerce in 1940. Inre- 1 
_ COrding and determining the general climate of the United States and its pos- eae 
_ sessions, the Weather Bureau ‘maintains weather, climatological, and hydro- ae 


logic services. A central repository | for its ‘data is provided at the National : 


4 Weather Records Center, North Carolina. Weather 
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a duly, 


| cooperates with federal agencies, its cooperative 

are limited by lack of authority to accept non-federal funds. neni 

- ‘The Weather Bureau has responsibility for forecasting weather and river 2) 

conditions, and for summarizing, analyzing, and interpreting related hydrolog- 

ic data for general use. It maintains a climatological and hydrologic r network i 
observation : stations the country. 


Program ‘Reappraisals 


‘og _ The Geological Survey i is 3 conducting a a comprehensive ve reappraisal « of the — 
effectiveness of its data collection programs, leading to possible adjustments - 
and improvements geared to present and future needs for hydrologic data, — 
_ The results of network reappraisals will be . published, probably in summary 
“form, 
_ Revisions are under way in the stream gaging network, which is et 
three principal elements: a primary network of long -term stations, 
secondary network of relatively short-term stations , and a third category 7 


"project operation or management stations. Increases are proposed in both 
_ primary and secondary stations to a total of about 10,000 active gaging sta- : 
Study continues on the stream sediment station network and the daily co: " 
= quality network. _ These analyses are expected to provide objective 7 “4 
and systematic approaches to the needs and improved criteria for establish - 
The Weather Bureau has its reappraisals bes begun 
to modify its ‘station objectives and locations, _ Two types of stations are in- — 
_ cluded in the climatological network: a a few long- -term “bench-mark” ” stations — 
~ selected for permanence, and a greater number of short-term stations. - Many 
_ observation stations are incorporated in both the climatological and the hydro- a 
logic networks. | The revised hydrologic network will be based on known ill 
of particular projects and investigations rather than « on area distribution. The ir. 
a stations are less permanent and can be shifted to meet changing conditions. aan 
Both networks are shown on available at offices of the 


The size and scope of programs for collecting hydrologic data are influ- 


a ‘enced by many factors which vary with | changing economic conditions. Among _ 
* the factors which may exert an influence on the magnitude of data collecting © ; 
are: availability of funds, usually = population | growth and | 


= stream pollution resulting from urbanization and i industrial ‘activity and 
such | new developments as weather modification and disposal of radioactive 


“a 


One of the most st important nt factors in maintaining an adequate data base is - 
af ‘the availability of sufficient funds. In general, more and better data could be _ 
we additional Legislators and administrators, however, 
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a, at both h federal and state levels, , frequently have not provided nae ap- } . 
x 3 = propriations for carrying forward a stable, comprehensive program of data ig 
collection, processing, , and publication. The results of operating with insuf- 
= atest funds are inadequate a areal coverage, discontinuous records, and de- i 
in quality of data. As measured in terms of what existing technical 
Mi abilities can provide, however, most of the publicly financed programs are 
nearly maximum possible results with the funds 


ca _ phenomenal, In 1924 the United States had about 9,800 domestic water supply i. 
me. systems, and in 1950 the number was close to 17, 000. (3) Industrial growth is 
even greater, and further expansion is predicted; industrial water use of 80 . 
billion gallons a day in 1950 is expected to be 215 billion by 1975. 4, _ Both public — 
and private enterprise in the water works industry need more data for planning 


new supplies and more efficient utilization of available water resources. , | 


q 


ion j es: Agricultural use of water increases as more land is brought under irriga-— 
ey “di tion and as supplemental water is applied to growing crops in humid areas. : 
a 4 The 1954 U. S. Census of Agriculture reports an increase of more than 30 
percent in the country’ s irrigated acreage between 1944 and 1954. Water use 
a | and control must be based on adequate hydrologic and agricultural ‘data. ae 
EP _ Additional hydrologic data are needed regarding water quality in most 
= a streams in order to continue to solve the pollution abatement problem. The © 

~ + federal government has a responsibility under the Water Pollution Control — 

of 1956, while the states have a primary responsibility for controlling Ae 
7 aa within their borders. . Both federal and state agencies depend on 

knowledge of river hydrology to measure the abatement ement programs 
is 


d to set standard lu nts. 

- * Even though adequate f funds may have been provided, the trained man- — ate: 
ol 44 ‘power to carry out the necessary investigations has not always been available. 
es oi _ This was particularly true during the war years, but it has continued beyond — 


recent war. Adequate staffing of the government scientific agencies, 


> Te 


q 


which are under Civil Service, has been hampered by salaries not coaapeiitiee, ia , 
with industry. Some segments of the engineering profession are advocating me. 
indicates need for accurate and verification of results. ‘The fa fact that 
4 reported results of “ rainmaking” often conflict with each other is an indica- ee) is 
tion that special data must be obtained in areas where weather modification _ i: 
is being attempted. _ Adequate information collected and analyzed in accordance — a 
a anticipation of the time when streams ¥ will be carrying new pollutants — 


Special c civil service rating in these categories as an improving the 
Factors 
‘The he increasing use of various forms of weather modification ae 
a with engineering standards, is the only reliable measure of the effectiveness r i 
¥ 4 and other substances not now present, early effort should be made to initiate 


collection programs which will establish base period of 
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= and the more critical the problem becomes the greater is the need for ne 
curate and extensive water data. Growing competition for limited water sup- 
3 p 


It is ¢ often difficult to separate the collection of nesteate data from the De, 
_ Processing, analysis, and interpretation. Observed data, for instance, may _ 
be of little value to the user unless made part of a more comprehensive study rf 
in which basic data have been related to other facts. The collector of basic 
fhe = Bas: has a responsibility to analyze and interpret. to a degree that the end 4 
q is useful in several waystoothers. = 
a The primary functions of most federal agencies which collect brdeolagie 
: data are not the collection of data but the operations which require the use or 
aia of data. Even the Weather Bureau and the Geological Survey, the 
most likely exceptions, were established for other purposes. Nevertheless, 
data are collected, interpreted, and published as an authorized function of the 
federal government in promoting the general welfare. It follows t that some bec, 
degree of data treatment also is a proper federal function. 
Some processing of basic data is necessary in order to put it into forms — 
usable by others. The broad general investigation and its direct results in 
terms of data are the primary responsibilities of the collecting agency. . The ¥ Bs 
adaptation of the data for specific uses is a responsibility of the user. 
a ‘It may he impracticable to define the extent to which Spey type ¢ of aia , 


processing is undesirable and excessive prossnsing is unwarranted. The ex- 

perience and background of the collector, which may be useful to the user, an 

are lost if only « ‘raw” data are presented. In contrast, too much | processing ae 


may infringe upon private waste public funds, and thus adversely 


Advances in the last decade have been in the field of process- 
ing than in the field of data collection. Standard methods and equipment for * 
CS various hydrologic factors were largely developed prior to 1940. 
New publication devices and media, however, are now making large volumes ul 
more readily available. Electronic equipment for 


Both the Geological Survey and the Weather Bureau are using sation . 
-_ electronic computing equipment for analysis of hydrologic data, Water con- 
trol and power operations | are e electrically computed by Bonneville 


- Opening up new uses of data and permitting new applications: J ee 
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-. Centralization of hydrologic data | programs has both ote and disad- 
vantages, An advantage to the collecting agency might be a disadvantage to “ae 
data user, and vice versa. Although conditions may vary excep- 


tions, the following are some of the advantages and disadvantages: 


y 

_ Records may be more uniform, consistent, and ‘reliable; costs may be 

lower because of less and higher efficiency; publication and dis- 


Early availability with a low degree of accuracy may be better than later _ 


and accuracy; of local interest are not adapta- 


devised to use the nation’s water resources. Appraisals of water policy, made ) 
= since 1950 by many official groups, mention the great need for adequate basic = ¢ 
hy data. In policy statements on conservation of natural resources adopted by are 
x national organizations of agriculture, business, and industry, numerous refer- 


1 ences are made to the adequacy of data. “These studies and reviews by groups a 
outside: the civil engineering profession indicate that the optimum amount of x 3 
basic data is not always available. . They recommend financial support to im- _ 

Pn trys data being collected today, to be adequate, must include not © a 


value. ‘The uncertainty of tomorrow’ s needs makes today’s a 
of opinion rather than a certainty. - However, the experience of the past gives 


_ Before reaching a decision as to the adequacy of certain. data, its relation- . 


"ship to need should be examined, standards and require-— 


that mere of data is not the real objective. 
= A critical examination of hydrologic data programs shows some over- be 


; onbhaaie on volume of data. The need of the user is not always for ‘more — 9 
data but rather for the right data which can be interpreted properly. — Those 
| concerned with this problem should ‘recognize the danger of demanding n more 
7 data asa substitute for professional ability to use properly what is available. ee 
Nevertheless, most authorities agree that both the quality of data and the 
Furthermore, attention should be 
given to with new needs and | dropping whatever is no 
The ciiane 7 six p incipal types of hydrologic aaa is viewed as follo 
the Task Group 
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Adequacy of Data Collection Programs 
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In the coverage of basic data are adequate, an d the 
_ information is readily available. Public and private records of precipitation 
Bt - throughout the country are well distributed and some are of rather long dura- = 
tion, _ The network of gages for collecting short-period rainfall intensity may 7 
be adequate. The greatest deficiencies are in the mountainous areas, in small j 
watersheds, and where unusual climatic and orographic conditions exist. “a 
The rainfall intensity -duration- -frequency reported in the Weather Bureau — 
papers has filleda long - -standing need. Progress also is being made 
in the study of areal distribution of precipitation from local intense storms. as if 
Further interpretation of hydrometeorologic data of this type will mae 
additional tool for determining runoff characteristics. 
_ Snow-water content observations are increasing gradually a and are riko 


erally adequate for operational purposes i in | most of the areas. Modest veal : 


ords and to increase the knowledge of snowfall and water content. prime oe) 


_ Except for evaporation data, the coverage is probably adequate from an 
_ engineering standpoint. Presently available data on pan evaporation are ade- 
quate except in areas of rugged topography. More basic study of the Lake ot 
_ Hefner type should be carried out to show water losses in relation to other 
climatic factors, so as to relate pan studies to open-water conditions. ‘ 
_--- The Weather Bureau reports that a reappraisal of evaporation data ona 


has increased and will continue to increase over a. period of 
. time as the demands for more adequate data expand. However, there will — 
“ never be e justification for measuring every stream at every significant change mye 
: in flow. ~ Any future expansion of the network of gages should be based on 
maximum possible use of correlations, employing data from key gaging sta- ae 
; tions to estimate (within close limits of error) t the runoff at other locations. 2 
Additional application of correlation techniques to present records ‘might - 
7 low some reduction in the number ot gages without jeopardizing accuracy. a 


of tidal influence, already so should be pi 
USGS stream gaging is now undergoing reappraisal to 
‘mine its atoquacy and the most suitable density to meet desired needs and 
mar _ standards. ‘This reappraisal is taking into consideration the views and needs 


_of other agencies and data users. Some expansion is taking place in analysis 


af ot runoff records which leads to or assists in later multiple use of data, 
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Ground water investigations have been expanded in recent years, how- g 
ever, and some of the deficiencies are disappearing. «. 


= 
- counties or - other political unit would improve the data r now being obtained. a r 

Exploratory work is needed to implement the discovery of new potential eee 

sources of ground water. As water uses approach or exceed safe yields, 


studies of limiting factors will become more detailed and data will becom =. 3 


_ more important. — _ Publication of information on the hydraulics of aquifers 
=". Although a network of federal observation wells is in n operation, the extent 7 
of its coverage is relatively small, Federal participation in data —— 


a the basic knowledge of ground water resources, ‘ 
Some ground water data are collected by investigators but Little 


appears in print. Only a few states require public filing of ‘well records, — ia 
water levels, logs, andpumpingdataa 


al in areas where it is increasingly important for 


The USGS is its sampling quality data in 


streams. The results may indicate the deficiencies and adjustments needed 7 7 


_ a The present sediment sampling program is somewhat less than compre- ae 
___ hensive, _ On the other hand, intensive sampling may be unnecessary after 
relationships between stream flow and sediment have been established. Some = 
sampling probably can be abandoned entirely. The study of the program now i. 
under way by the USGS should show were readjustments can be made to give 
PS Data on the accumulation of sediment in reservoirs are generally adequate i 
7 and should be continued. Small reservoirs are not as well investigated me A: q 
large ones, however. No one federal agency needs to be given primary re- 
G sponsibility for this work if there is adequate coordination among agencies, a 
fa] 


a of all the basic Dediantonte data collected in the United States, only a por- os 


_ is prepared for publication in reports, The regular publication media of 
7 the | Weather Bureau and Geological Survey are 1 more complete | than they ywere ya 
a 20 years ago. In general, they are meeting the needs of the engineering gal - Ky, i 
grams - especially those involving large river basins and major drainage 
areas, on Water Resources has compiled 


24 
4 
| 
4 
i nh chemical quality and temperature Ol water 1 

more often due to lack of uni and standardization of field practices wy 
of investigation. Better co 
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useful inventories of data on water quality a and sediment load. 
‘In recent years the demand has increased for hydrologic data i in nthe > smal~ 


7, :: ay issued which help to meet some of the demand for this type of data. ree 
Many private collectors of hydrologic da data have no plans for publication of 
_ their basic data. a There is no established procedure to have ‘such data in- — 
& cluded in regular or special publications of federal or state agencies. Data - 

; thus collected for a special purpose may not be suited | for publication in in jour- 
nals of scientific societies « or in special bulletins of private or public organi- A 

_ gations. In such form, also, the basic data might not be sufficiently complete . 


to afford a check on conclusions which may be drawn. Much unpublished data — z 


a filed upon . completion of the study, and is either lost or destroyed, 
_ _The time lag in the publication of precipitation data has been greatly re- 
duced in the past f few years, , now being generally le less than two months for ; 
hourly and daily precipitation. Stream flow data appears s in print a about 24 a 
-months after the end of the water year, the delay being due to the time re- 
a for analysis and review and to a backlog of computation work in some 
field offices. New p procedures being considered for stream flow records 
( fa to shorten the lag time for publication by several months. " Decreas- 
ye ing the present lag time for precipitation data does not appear justified. pret he 
Descriptions of the regular and special publications of all the federal weil 
+ agencies dealing with hydrologic data are given in the previously noted report 
*Principal Federal Sources of Hydrologic Data.” 
‘The periodic and serial publications of most agencies are improving, but 
the need seems to be for different types of reports. There is need for all = i 
special summaries, such as those issued by the Geological Survey an and the the iad 
Weather Bur ¥ Some should be issued more 
4 8 It has been recommended that some means be p provided for publishing ng data 
= _ collected by those who have no publications program. Examples are the vast | 
B ys amounts 0 of f sediment data collected by th the Corps: of Engineers and the infor ae 
a mation on water quality collected by numerous agencies. The volume of un- | a 
i published hydrologic data constitutes one of the profession’ svaluable 
= resources. How to make it more readily available and useful is a challenging — 
problem. The designation of selected agencies as official repositories 
as The public agencies now devote large portions of their money and man- 


power to preparing data in forms most useful to practicing engineers. They 
- constantly seek to produce better and more useful publications. If what — 


prepare is the fault may be only partly theirs. the need is for 


The new high- -speed methods of | machine are of proven ave in 
“tet 


cards or data available in in this fields been 


| 
a ___ develop sound engineering plans for watershed protection and flood prevention —_ Bi. 
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or the data collection agency to place the data on machine punched 


= 


"HYDROLOGIC DATA 


opened for ready processing of analyzed or interpretive data by the collection 


agency. The cards or tapes can be used and reused for this purposes 
aa Advantage may well be taken of high-speed machine ‘methods to increase — 
the he amount of analyzed or i interpretive data to 0 be i issued | by the collection 
agencies. " The limit of such analysis and interpretation is difficult to —. 
lish. However, the data collection agencies should not unduly emphasize the 
and interpretation at the expense of basic d data collection. The fol- 


Collection of basic data should take priority over analysis and interpre-— 
tation. Analysis and interpretation can take place at any time, but 
a ; most cases hydrologic data can only be collected at the time of the physi- 
ate cal occurrence. Recognizing that a practical limitation will always be _ 
put on the amount of funds available, analysis and interpretation 
not be ‘overemphasized so that the basic data collection program is 


lowed to fall below a generally acceptable 


: ne Analyzed or interpretive data should be processed if there is reasonable = 
Pie expectancy of subsequent multiple use, but should not be processed _ 
to serve a single purpose. 


3) In most cases, the collection ag agency is best qualified to make interpre- 
i _ tations. The basic data used in the interpretive analysis should ~! 
available to the user for his own interpretation. 
Within these general limitations the collection agencies should 


7 ‘The high- pre methods of machine computations have additional implica- 
7 tions in the data collection program, The use of electronic computers for — 7 
ae analytical and interpretive purposes can be expected to increase, and as such : 
use increases, the ability of the user to consider greater amounts of data bie 
a will correspondingly increase, The result could well be a demand for greater 
G4 amounts of data, whereas in the past such demands may not have been exerted 


because of the of processing voluminous records ; by manual 


4 is frequently Decision on whether or "not to use microfilm 


_ must be made by the data collection agency on the basis of its own knowledge — 
e of the probable future use for the original data. As one of many factorsin es 
_ the internal efficiency of the data collection agencies, it would appear to have “i 
A ‘ry on the adequacy of eee data only to the ‘oment that efficiency 7 ‘ 


wae q 
= 
= 
rainfall intensity probability, streamflow frequency, summaries of 
selected occurrences such as maximum streamflow or precipitation, and 
| _ The data available on cards or tapes have also been used in other ways. any oe 
. _ Copies of cards or tapes can be made for the user who can then program, on © ' 
aS a machines available in his own locality, the analytical processing he requires, — Bee 
™ _ Fuller use can be made of cards and tapes as engineers become aware that 7% aa 
4 
P 
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ba: new Types s of D Data 
ost authorities : agree that water 

and more places in the years ahead. All trends of water use point to intensi- 
fied competition for water and ways to use it. Few engineers 100 years ago i 
appreciated the need for maintaining records of stream flow, but none doubts 
the value of such records today. Twenty years ago no one thought of testing 
our water courses for radioactivity or to determine minute amounts of trace 


sg elements. Such tests are only now becoming important to engineers and 

g ha received on the e earth, most o of our hydrologic-« engineering problems —- 

a ae better use of what we have, using it wisely and treating it well. For example, — 
evaporation generally is a waste of valuable water, and if evaporation can be 
p reduced the effective water supply will be increased. ? Encouraging results a 


_ have been obtained from experiments with the use of monomolecular films as 
"inhibitors of evaporation from water surfaces. 
other fields as well, ‘such as” sewage reclamation, | new opportunity | 
be found to develop good water supplies when and where needed. _ The recla- 


be learned of the potentialities in the United States. So far reclamation has 

a _— been developed principally in areas of critical water supplies where it is ; 
feasible to process a large volume of sewage f for industrial or irrigation use. 

Hydrologic data are needed on both quantity and quality of polluted waters, q 

including not only domestic sewage but also saline (natural) waters which = 

might be reclaimed, _ New: data may be needed v when the — =a 


The use of radioactive materials in manufacturing, medicine, a and ithe 

a production of electric energy, as well as nuclear reaction f for defense pur- 
poses, will increase water pollution unless great care is exercised. ‘The ae 
lector of hydrologic data is one of the watchers in this new field. New data _— 


fects of watershed treatment measures. on the physiographic, 
climatic, and other conditions, ‘such measures m may have various effects, | 
_ hydrologically and economically. © More data are required to assess such ef- 
Ba fects for di different parts of the country. . Itis noted that a joint study to > develop 
iS. ‘3 procedures for making ng determinations of hydrologic effects is being carried an 
out by the Agricultural Research Service, Soil Conservation Service, and 
Loa Artificial nucleation of clouds to increase precipitation has been the si 
ee of research and actual operation for more than a decade. Arecent report 
H ‘ad indicates that at one time more than 20 percent of the United States was being _ 
subjected to rainmaking efforts. Commercial firms claim considerable suc 


cess, and beneficiaries willingly pay for their services. In contrast, however, 
(oe World Meteorological Organization reported in August 1955 that a net in- 
crease of precipitation has not been demonstrated beyond reasonable doubt in 
f ‘any seeding operations yet described in technical literature. Yet, it is in- = 
i correct to state categorically that rainmaking is impossible. Much orderly y 
Ae research will be required to establish the potentialities of weather modifica- — 
tion, ‘Statistical treatment of hydrologic data is the only way in which the re- \ 


sults | of operations can be assessed 
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APP, 
HYDROLOGIC DATA 


data, long- ‘range planning ‘should to determine new needs and 
types of data. It seems obvious, for example, that more complete andac- 7 ; 
curate data is is needed where prospects of water problems are | most iminent. we. 
‘The civil engineer should exercise the imaginative thinking necessary to if 
provide a capacity to collect and interpret hydrologic data as" will be equal | 
the sum expanding demands. iy 


Ory 
The following have been drawn the by this 


. The increasing urbanization and industrialization of the United States is 
ec ™ continually greater amounts and better quality of hydrologic data. a 
Considerable improvement has taken place in recent ‘years, particularly with- 
_ in some of the federal agencies, but further progress is necessary. In gen- 
eral, adequate hydrologic data are being collected to meet most engineering 
“needs. AS deficiencies in basic data are overcome, , however, more emphasis 
- must be put on analysis ‘and | interpretation in order to increase the he availability ; 


2. Deficiencies continue to exist in certain types of data: precipitation 3 

and runoff data in small drainage areas apart from agricultural and industrial Fs 
_ developments; ground water data in undeveloped and newly developed locali-— 
ties; and water quality | data in a number of scattered regions. Pe, HED stor 4 


_ 3. Publication of many types of hydrologic data is improving and the time — 
ia are decreasing but not in proportion to the growing needs. New ki kinds of 
_ reports, such as periodic summaries, are making some of the widely | used till: 
data more readily available. New machines and processes are capable of re- 
ducing the time (but not necessarily the cost) required to produce usable A > sata 
records of basic data, thus leading to greater epee by professional 
both in processing and in using the data. 


4. Ifa suitable means could be litte for processing and publishing mis a 
7 ‘ eiimiene data from public and private files, valuable records now lying use- 
= be made available to the the engineering 
“Inadequate federal and state appropriations have held back some ‘worth- 
- data collection programs, Progress is hampered by lack of qualified _ 


professional personnel who are not attracted by government salaries which — / 


existing federal networks of gages Measuring hydrologic data are 


= about improvements to eon pace with changing needs and conditions. 


failure | to begin re records needed for the future. 
cooperative programs between state and federal are 
yielding satisfactory data on streamflow and ground water. The extension ‘call 


——- programs to other types of data might help to retain state - 


agencies. 
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8. The of federal and state activities in data collection 
regional interagency committees is desirable. Similar coordination tog 


federal activities among various agencies should be encouraged. igh 


the engineering « competence r required to interpret and apply the data | for bene- _ 


recommendations 
Federal- state cooperative programs of basic 


should be extended. In addition to ‘strengthening the cooperative ‘programs 


for stream flow and ground water, similar | programs should be established bai : 

for other types of hydrologic data, such as hydrometeorological, as recom- 

mended by the President’s Advisory Committee on Wi Water Resources Policy. 
Sach ‘Programs should be encouraged : at the state level. 


2As satisfactory federal hydrologic data program can be obtained elements. 


in a single federal agency. 


Sound conservation and development of the nation’s water resources 


quires more adequate hydrologic data which can best be derived by applying — a 


the technical abilities and the professional integrity of civil engineers and — i 
3 others with similar qualifications. Therefore, competent persons should be ay 
in charge of designing hydrologic data programs and supervising the collec- Ea 


tion, processing, and scientific use of hydrologic data. 
g a 4. Having both a professional and a civic responsibility with respect a 
that data n no longer needed for private would be availa- 


_ qi hi hydrologic data, the private civil engineer should maintain adequate records 


7 a me conclusion, the Task Group appreciates the opportunity to serve the 


‘Hydrology Committee and the Society. The Chairmen sincerely ~ a 
a oe the assistance and cooperation of the members of the Task Group, the mem- 


Resp 
Robert R. ‘Balmer Charles C. “McDonald 2s, 
©. O. Clark Edward J. 
Clayton H. Hardison Walter G. Schulz 
Lloyd L. Harrold William W. Wheeler 
= Herbert E. Hudson, de. Donald W. Tan Tuyl, 


‘Task Group on Hydrologic Da Data of the Committee o on 
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“HYDROLOGIC DATA 
REFERENCES 


sponsored by the Conservation | Foundation, is expected to to ‘published 

per 

1958 under the title Data for Decision. 


E 2. _ Federal Inter - ~Agency Committee « on Water Resources, Principal Federal — 


P= I — STATEMENT BY U. S. WEATHER BUREAU ~ 
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four years ago and is now completed, The reappraisal had two purposes. One =f 
was to determine the required density and distribution of precipitation and 
temperature observations for adequate definition of the climate. ‘This deter- 
mination has been followed by authorization of a representative gridofob- | 
gradient and 
< 
necessary stations are now being discontinued, and new stations S established 
in accordance with the plan and available funds. 4 
a A limited number (25-50) of stations with long records aaa nei 
_ pects of continuity will be designated as bench- mark stations, and ~" > records, 


a a The next step ii in the climatological network (still i in 1 the pinnnles stage) is Bs 


to utilize the double category climatological network plan. Such a network _ 
_ would consist of one group of stations which we expect to keep unchanged al- 


> most indefinitely, and another larger group, each of which would be moved ito 
©. every 5 or 10 years to provide data from a greater number of locations, alee 
a from these short-term stations could then be correlated with data from near- 
by long-term stations to provide an estimate of the long-term climate for the 
_ The other purpose of the reappraisal was to establish - known Pe: 


the needs for both and operation in the field of applied 
large part of this requirement was based on the 1949 recommendations of the 
_ Federal Inter-Agency River Basin Committee. ‘This hydrologic network in- $ 
cludes portions of the climatic network but not follow general rules 
distribution and density of observations. It is based, rather on the 
needs of particular projects and investigations. — It is, therefore, flexible and 
to change; and its implementation depends on project authorization 
availability of funds for these projects. 
iB — networks are shown non state maps s used as working tools in our sub eat 
station ‘management pro program. _ Publication is not intented but they can be 


inspected State Offices and copies made available if required. 
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determine the climate of its possessions, , at the 
. same time making available the basic observational data by archive, punch ~ 
card, microfilm and publication, In addition, the Bureau has extensive fore- 
7 1 casting responsibilities, both for weather and river conditions, for which ad- 
oh ditional observations are required and these, too, are incorporated in its me 
basic data program. For the more specific studies, both climatological 
hydrologic, required by c other agencies and specialized interests, and 
- operational use, the Bureau provides its facilities, its advice, and its co- 
fame whenever possible. For all these data, and for other data made 
available to it, the Bureau provides the national repository at its National | 
e.. eather Records Center at Asheville, North Carolina, = = 
f ab The extent of summary, analysis, and interpretation of these basic data is 
to not spelled out, except that the Bureau’s s self- -supported limited program io * 
along those lines is naturally for general rather than project use, unless the ~ is 
project is the improvement of the Bureau’s forecast procedures or climatic 
; presentations, _ However, as in the observational program, the Bureau also 
; ie provides its facilities, advice, and cooperation to meet other agency or inter- 


agency needs in ‘Summary, analysis, or interpretation. numerous 


Papers Nos, 2, 4, 10, 13, 15-17, 23-26, 28 and 29, and Weather Bureau Re- > 
_ search Papers Nos. 34 and 38 are cited as examples of what | can be done re | 
es. funds are available; in a majority of these cases the | costs of these 


cies. 


vancement of basic data programs, both and analytical; and 
has been facilitated by inter-agency transfer of funds. same 
facility is not available to the Weather Bureau for cooperation with State or 
other non-Federal organizations, although the Presidential Advisory Com- | 
mittee on Water Resources and Policy has recommended extension of — ae 
_Federal-State fund- -matching and similar authority to all Federal basic data ee 
agencies. Certainly extension of this practice to the meteorological and hy- ad ai 
& meteorological field would stimulate State interest and permit intensifi- Pt 
cation of networks and special investigations considered to be of insufficient — 
a - interest to justify Federal expenditures. The purpose, as in Federal 
cooperation, would be to realize | the technical and economic 
form required by the design specialists of me agencies responsible for de- 
operation of water utilization, ‘conservation and management 
ie 


Data programs under a single agency ene 0G 
; ee Centralization of basic data programs within a single 2 agency should provide _ 


ee the usual advantages of more uniform standards, closer coordination of re- 
lated programs, , and elimination of duplication. However, coordination and 
cooperation among the agencies concerned are now very good, 
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21 
Witness, for increase in ‘multiple -purpose 
‘serving several agencies; transfer of funds from agencies like the — of Bx! } 
f ‘ Engineers, , Bureau of Reclamation, and Soil Conservation Service to the 
f _ Weather Bureau for the establishment of meteorological stations and for 
special hydrometeorological studies; similar transfers to Geological Survey; 
and cooperation between Weather Bureau and Geological Survey in prepara- ; 
1 of hydrologic atlases and flood t 
ion of hydrologic atlases and floo repor s. 
= ae is thus a natural, healthy growth toward realizing the advantages of 


and aiiaianel ‘supervision of diverse aspects s of the basic data. program. 
_ The meteorological and climatological data of the Weather Bureau, for ex- p= @ 
i ample, are collected fora a great variety of uses, , of which only some are hy- 


drologic. _ There is no more reason why they : should be collected by a hydro- a 
; logic data agency than that they should be collected by an n agricultural, aero- 
% nautic, ‘health, manufacturing, or marketing data agency. . They are essentially 
meteorological data and are likely to be most efficiently and 
collected and processed by a basically meteorological agency. Even the 
_ Bureau’s most t strictly hydrologic data, such as river Stages, are so intimate- 

_ ly tied in with ‘meteorological observations and forecasts for river forecasting “ 

_ purposes that they are most logically collected within the same channels and 

under the same auspices. . The single ~agency concept, it seems to us, arises 

rases 
from an n over -simplification of the problem. A more realistic ‘approach would — 

be to consider a regrouping of elements of the basic data programs of on 
various agencies, in order to bring together elements ‘most ‘Closely in- 


terrelated within the appropriate agenc eae, 


This could well be increased and is another heading. 


factor to guard against in such cooperation is the year- by-year uncertainty of 
financial support that epee jeopardize the wi of records and make 


STATEMENT BY U. S. GEOLOGICAL SURVEY 


Reappraisal of networks 


= In 1955 the Geological Survey began a comprehensive study of the network 
_of river measurement stations, chemical quality sampling stations, and omit 
ment measurement stations in ‘the continental United States, Alaska and 

Hawaii. These networks of stations together with the ground-water observa- 
- tions wells constitute the Survey’s continuing program of systematic meas- — . 4 
urements and observations of ‘te occurrence, movement and quality of the a 
_ water resources of the nation. The study was designed to determine the ef- “— | 
 fectiveness of the network in defining the areal variations in the factors being 
measured and to serve as a guide in making adjustments < and improvements a 
to help meet current and anticipated needs for data. 

_ The network analyses now being carried on are merely or one step in ria con- — 
of the of the programs of collection of | basic hydro- 


i 
2 
| 
4 
Ae _ many respects an integral and inseparable part of, the total investigational oo 7 


: basic measurements and observations. This is particularly applicable to in- 
i vestigations of the character and availability of ground water. The ef ne 
nation of the geologic, , hydrologic, and hydraulic characteristics of aquifers = 4 
a complex problem requiring close control of investigations and on-the-job 2 
analyses and interpretations of data. Networks of measurement stations © ail omit 
; therefore are only a part of the total objective of providing adequate knowl- 7 
edge of the water resources, results of these network appraisals will 
be published in appropriate forms as promptly as practicable. The great ese 
amount Of detailed work involved indicates a summary form of { publication, 
‘The stream gaging network is composed of three principal components. be 
primary network of long-term gaging stations, well distributed geo- 
Bilis graphically and operated indefinitely, provides an index of the long- range 
_ fluctuations in streamflow for any region, A secondary network of relatively 
short-term stations defines the areal variations and characteristics of runoff; 
data from stations in this secondary network and other ‘data collected on a | a 
short-term basis may be correlated with data from the permanent type pri- 
mary network, thereby achieving perspective in relation to long-term climatic 
he and hydrologic fluctuations. Stream measurement stations that provide data 
specifically needed for project operation or water management purposes 
comprise the third category. The primary and secondary networks are sup- 
_ plemented by measurements of flood and drought flows and other oe 
tions as required. Some of the technical aspects of this network analysis have _ 
‘The study of the stream gaging network is essentially | complete, and re- 
i a - ports for various states and regions are under preparation. While the study _ 
a c cannot be expected to provide a complete measure of the deficiencies in data 
_ and network requirements, it does contribute substantially to the appraisal of 
current needs for data and provides a guide for future network expansion. It a 
_ indicates, for instance, that the existing primary network of 2,600 stations 
_ should be selectively. expanded to about 3,300 stations, and that nearly 1,300 i 
existing stations can be relocated after a period of operation of 5 to 10 years | } 
duration. The study also indicates that about 1,150 previously discontinued _ - 
> iy stations had accumulated sufficient records to permit adequate correlation PF 
base stations. Also, 1,250 existing stations operated for water manage - 
‘iz 4 ment purposes also qualify for the secondary hydrologic network. It was con- © 
. cluded that about 2,500 additional short-term secondary stations should be oe 


program which includes the processing, analysis and interpretation of the 


: { established promptly at other selected locations for a well-balanced and cur- 
rently adequate stream gaging program. A total of about 10,000 active gaging 


stations is indicated as the desirable network size to meet present require - 


ments. _At the present Gme. 6,640 stream flow measurement stations i. 


Daily ‘records of sediment movement in streams are obtained at approxi-— 
ae 150 locations. Because the coverage is relatively sparse and large 


“Thest study is not yet complete, but activity to date indicates that it will reddk a 
in improved criteria for selecting (a) the optimum composition of a sediment iis 
with respect to ‘measurements, analyses, interpretation, and re- 
_ search at a selected level of activity; (b) the distribution of observations with 
respect to time and place; and (c) the characteristics and properties of s¢ sedi- 


‘ ment and its watershed environment that should be observed to maximize 
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Er As of June 30, 1957, there were 300 daily chemical quality stations me - 


gg better understanding of the physical processes in the interrelation of the a xi 
fe various factors influencing the occurrence and movement of water facilitates £3 


Ee _ investigations, and oy criteria for the establishment of stations com- 


utility c of the of the 1 is underway and 
complementary studies will follow. At the present time the broadly stated ui 


mately the s same that for the sediment being limited 
scope because of the relative sparseness of the coverage. The analysis is pot 
ie expected to provide an objective, systematic approach to chemical quality — i 


‘major in his approach, ver 


The broad objective of the ‘Survey’ s program with respect 


_ water resources is the determination and description of the quantity ee 
A quality of the surface and underground waters of the Nation. ‘This broad ob- . 
2 jective encompasses the appraisal of the total oer resources; the variations 


ots water through the permeable portion of the earth’s surface; ant the ‘mag- 2 
nitude and frequency of floods and low stream flow. OF 
_ For many years the water resources problems of the country were a: 5 
‘3 that primary emphasis in the Geological Survey’s program was placed upon a 
os the collection of basic water resources data. The increasing demand for _ * 

‘Solutions of increasingly complex water problems now necessitates greater ps! 
attention to analytical, interpretative, and research phases of water resources o 
investigations. The Nation is at a stage of development wherein a concerted — te 
oa attack is urgently needed upon all phases of water resources problems. The pid ; 
_ Survey recently began to increase the amount of work being done on those 
_ phases. Increased and intensified effort is also being directed toward re- oat 
search in the basic physical processes involved in the occurrence and move- 
ment of water in order to understand more fully the various facets of eo 
phenomena which we are observing and interpreting. Experimen- 
7” tation and additional research directed toward the development of improved — = 4 
techniques for investigations and analyses have been carried on and are now : 
a in progress by the Geological Survey. The investigations include study of 
movement of water through unsaturated media, evaporation from free water | 


‘surfaces, evapotranspiration, the effects of land use, sediment transport, 
channel geometry and river morphology, erosion, infiltration, and recharge. 


more effective collection and interpretation of basic data. Consequently im- 4 
j a knowledge is available | for application to the planning and design of © ies 
Adequacy o of | basic enabling legislation 


F an The Act of March 3, 1879, "creating the. Geological Survey, , did not define 
as responsibilities in the weer ‘resources field. 
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- ated funds for water resources investigations, In the 1894 statute, Congress» 

made funds available to the Geological Survey for “gauging the streams and _ 
determining the water supply of the United States, including the oe 


3 appropriation acts, ‘1888, ‘have authorized the use > of 


‘a of under-ground currents and artesian wells in arid and semi-arid sections.” 
Successive appropriation acts have continued or further define the authoriza 
Le tion for the water resource functions of the Survey. The current programs — 
_ and any presently contemplated future recommendations are fully authorized an 


’ - current trend is for increased participation by states in cooperative water se 
resources investigations. i The Survey will continue to acknowledge this par-. 
ticipation to the greatest extent possible within its available financial and 
_ manpower resources and in accord with the desirable intention of maintain- 
ing a balanced over-all program of mapping, investigations and research. eT. 
Because the programs are geared substantially to water 


4 crease in water development and use, water problems are 
oe localized and more regional in nature, requiring collective action on the part i 
4 of the States and increasing the need for studies of the water resources ona 
broad geographic basis. This condition has —_ the grember need for 
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a = _ For more than 50 years the mutual interests of the State and Federal 7 Bix 
— governments in water resource problems have been implemented in the 
Geological Survey by cooperative investigations. These cooperative activities She 
been majo part of the work of the Survey, and through them far more 
ha sic data have heen obtained than would have heen acquired othe ise =§The 
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GROUND WATER SUFFOLK L.L, NEW YORK? 
John F. John M. 
a 


ot 


Ground water pumped from the aquifers underlying Suffolk ng 


‘a 


Island, N. Y, for all uses: amounted to more than 24 billion gallons in n 1955. 


a clude: extensive sea water encroachment, marked as of storage, in- 
ae contamination, and rises in water temperatures. bani 


count Ty one the growing in the 
"States, has an area of 920 ‘square miles — — about three- -quarters | as large as 
| Rhode Island. As can be seen in figure 1, it occupies the eastern two-thirds _ 
_ of Long Island, with ‘Nassau County, another 1 rapidly growing county, on the a Mr 

* west. Sea water bounds the county on three sides, dividing it geographically - 

A into three areas: the main body and two peninsulas, termed the North and = 
a South Forks. "These fo forks ‘ks originate 2 at Riverhead and extend eastward for 

- about 27 and 40 miles, respectively. ‘5 Although creating a Sevosniie environ- 


hy County’s ground- supply. Suffolk County’ s ‘rapid industrial ‘development 
has been, in large measure, due to the abundance and quality of this supply | oo 
and as continued industrial growth at an increasing rate is a strong possibili- 
i. ty, investigations of the hydrology and geology of Suffolk County in relation to 
sea-water encroachment are of utmost importance. 
~ The purpose of this paper is to outline briefly the ground-water pieibene Q 
_ and geology of Suffolk County, to describe the saigsian involved in oem: 


a. Publication authorized by the Director o of the U.S. Geological ae 
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Pca ok water fins the pore spaces | of the > unconsolidated clays, sands, and 


The lithology of these aquifers is s briefly described telow. More 


coverage is made by de Laguna and Perlmutter (1949), 
i om The Raritan formation, deepest of the unconsolidated deposits on long 
Island, rests on a basement of crystalline bedrock. Raritan is composed 


of two members - an upper clay member and lower Lloyd sand member. 


i about 80 feet per ‘mile ina _ southeasterly direction. The clay member beneath 
' most of Long Island usually consists of beds of silty and solid clay with oc- ce 
- casional sandy layers. The Lloyd sand member, on the other hand, is a pre- 
& dominantly coarse sand and gravel intercalated ‘at intervals with thin layers 
- of silt and clay. As the overlying clay member has a very low permeability, — 
the water in the » Lloyd sand member is for the most part confined under b, 
4 artesian pressure. The Magothy (7?) formation is considerably thicker than _ 
. the underlying Raritan formation in most places on Long Island, It consists 
of layers of fine sand, silt, and clay interbedded with several zones of coarse — 28): 
- sand and gravel. Pre-glacial erosion of the Magothy (?) formation has caused at 
5 a considerable relief to be developed in the upper surface. Consequently, the | a ot 
- depth at which this surface is first penetrated by wells is quite variable. The 
4 - water in the Magothy (?) formation is generally under artesian pressure, but | 
locally it is unconfined. Glacial deposits of late Pleistocene age mantle the 
older formations in practically all places on Long Island. They form the 
=. present land surface and surficial deposits of Suffolk County. Two different _ 
types of material are recognized from the cuttings of wells drilled into the 
i _ deposits - till and glacial -outwash, Till is a heterogeneous mixture of sand, ps 
= clay and boulders deposited from glacial ice. It is generally poorly er 
ms and the presence of clay, in some places, causes this material to have red 
low permeability. In other places, however, the till is essentially free of clay 4 
_ and quite sandy and locally may have high permeability. Outwash, onthe } 


other hand, was 8 deposited from glacial melt water and is moderately to well 

sorted. It is largely composed of sand with some e gravel and has a relatively 
high permeability. Stratification, however, and the presence of lenses of wh 
and silt, , interbedded with t the sand and gravel usually (Conses the ho horizontal — 


permeability to exceed the vertical permeability. 

_ _ The water table of Suffolk County, the upper limit of the entire thickness of : 
_ saturated material, , ranged in altitude from about ‘10 feet above sea level in 

a the inland, central portion of Suffolk County to sea level near the shores. The e 

total amount of water stored below the water table amounts to many billions ie 
_§ gallons and depends chiefly on the volume | of saturated material and the Ba! a 


porosity. The aquifers are constantly receiving replenishment and are con-— 
Ss discharging. Variation in either replenishment, discharge, or both 
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_ shown in figure 2. Three aquifers have been recognized. From bedrock up 
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mean solely from precipitation, which averages about 43 inches per year. 
Only part of the precipitation reaches the water -bearing sands, for sizable 
- amounts are lost by evaporation. Losses by overland runoff to streams in A” a 
undeveloped areas are small, the rapid of precipitation 


- cent of the precipitation reaches the water table or an average of about 350 ci 


‘Natural discharge from Suffolk County’s 3 ground- -water reservoir takes 


the water table. According to furnished by R. M. 

_ Sawyer of the U. S. Geological Survey the average discharge by cosag- coal 

< is about one- ~quarter of the average ground-water recharge or more than - 

80,000 million gallons annually. is about 3- 1/2 times the quantity 
ground water pumped for all purposes in 1955. Likewise, an unknown but 

. "probably greater, quantity of ground water moving through the subsurface } 
; sands, gravels, and clays is discharged each year directly into the ocean and 
into streams below the gaging stations. A third water loss, also unknown, but 
"possibly sizable in amount, occurs through evaporation from surface-water 
bodies, such as streams and ponds, and evaporation and transpiration from 
the marshy and low-lying areas that fringe the shoreline of Suffolk County. — 
‘These are the natural elements of ground-water conditions in Suffolk 


«: County. _ Ground- -water development can markedly alter the natural pattern. 
- Some of | the changes have a beneficial result, others ‘create problems. ; 


ial Water utilization, in general, requires water adequate in n quality as well as 
“ in quantity. Substantial amounts of dissolved solids in the water can cause a - 
number of problems. Scale in boiler tubes and piping, resulting from the de 
pos sition of the carbonates of iron, , calcium, and magnesium, impairs heat a 
_ transfer and increases head losses. Encrustation of well screens by these or 
other compounds reduce both the yield of producing wells and the ability of — 
echarge wells to take water. | Sizable. concentrations of iron and manganese 
an stain laundry and water fixtures. Excessive concentrations of nitrate, _ 
ie: sulfates, and chlorides can produce organic disorders. Water with a pH 1 less 
than about 6.0, especially in the absence of carbonates, has a marked corro- = 
sive effect on ‘metals, causing equipment failure and deterioration. Accepta- _ 
a ble standards for water quality differ depending on the particular use desired, — 
. Ih general, a water supply that conforms to the chemical and bacterial stand- 
_ ards for drinking water is suitable for most purposes. 
a Figure 3 shows the locations of 1b water - ~supply plants in Suffolk County. i 


_—— 6, 000 million gallons of ground water for public supply, laboratory, and 
institutional usage (see footnote - table 1). According to estimates of the = 
New York Water Power and Control Commission this pumpage is about one- 
2 quarter of the total ground-water pumped in Suffolk County in 1955. “ee 
‘Table 1 shows the type of water treatment employed at these wells to = 
water of satisfactory quality to the consumer. The data were obtained 
from the Chief Engineers or the of the respective No 
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treatment whatsoever is necessary at six ix installations. _Corrosion control is 
“necessary at sixteen well installations and is the only treatment widely used. 


rather than a remedial measure required by the State Health Department. Bs. 
Treatment for iron, turbidity, and odor is needed only in isolated cases. Be- 
cause of the low concentration of carbonates and of magnesi- 


conditioning must have a relatively low temperature. Here again, Suffolk — 
a County is fortunate, for in general, ground- -water temperatures range fr from ; 
50°F to slightly more than 64°F 


Le 


In 19 
a In 1955, according to estimates of the New York State Water Power ho 


ation and the Lloyd sand 
member of the ‘Raritan formation is relatively 
Ve 


Methods of 


‘bearing sands underlying Suffolk County. For the most part, these range in in 

depth from a few feet to about 800 feet and have steel casings ‘ranging in i 

_ diameter from 1-1/4 inches to about 12 inches. Well screens which are a — 
necessary part of the well construction range in length from 2-foot ‘well-point 
“screens on domestic wells to 75-foot (or more) screens on major public- cae. 
Supply wells. An important part of the well drilling operation is well develop- te 
“ment, which involves the removal of the finer-grained particles from around 


hs _ the screen. Were it not for well development, many wells en in uncon - 
solidated material could be pumped only at very lowrates. 

_ Where the depth to » water is less than about 8 feet, artificial ponds are oc- 
-casionally t used as a source of irrigational water. ‘The dimensions of the 

ponds average 40 feet and 6 feet deep but differ for 

-well are the most commonly used type in 

. County. Centrifugal pumps, although used to some extent, are limited by 

re = their inability to raise water continuously from depths greater | than | about - el 
- feet. _ A modification of the deep-well turbine pump that has recently come into 


q ag popular use is the submersible pump. All the working parts of this type pump, 


- including the electric motor, are contained inside the well, being submerged 
below the water level. ‘Where small -capacity pi pumps are e required, ‘such as 
“ for household use, and where the depth to water is —— suction lift, jet 
Pumps for irrigation wells are e usually driven by or diesel en- 
ot for public-supply, industrial, and domestic wells are usually : 
_ driven by electric motors for smoothness and efficiency of operation. Most — = 
—_ ‘supply and large- capacity industrial wells are equipped with a gasoline 


or a diesel engine for stand- by use in. case of power failure. ea "Tac 
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Publics -Supply ‘Usage 


Withdrawals for pu public su supplies, which ‘accounted for three-tenths 


t the total pumpage in 1955, are made at 27 municipally-owned and about a E 


, 


privately-owned water- -supply plants. ‘The largest system is that ofthe 
_ Suffolk County Water Authority which ‘supplied about three-fifths of the ground 
_ water used for public supply. Water pumped for public supply is used for a 
variety of domestic, industrial, and agricultural uses. In eastern Suffolk — a a 
_ County the uses are predominantly domestic and irrigation; in western Suffolk 
County the uses are domestic and industrial. 


velopment of a community. Ins some areas s development may be due to indus - 
tr y; in others, to agriculture. Community growth, both in size and wealth, 
results in increased water use owing to the large number of water-using de- 
vices. For purposes of this paper, figures for daily use per capita in Suffolk — if 
x County were obtained by dividing the total plant pumpage in an area by the “S 
number of users served by that plant. New York City, a highly commercial- _— 


4 ized and industrial community, has a daily use of about 130 gallons per capita. 
A rural community with few water - peat devices may have a daily use as low 

A _ Listed in table 2 are the 1955 daily use per capita for some water - -supply = 

_ plants in Suffolk County. The locations of these plants are given in figure 3. ay 

_ Figures shown in table 2 are based on data obtained from the Chief Engineers | y 

or the Superintendents of the plants involved. The highest daily use of 126 ‘< a 

gallons per capita was in the Incorporated Village of Greenport, a developed, 

a sewered community. ‘Here, more than one-tenth of the pumpage is used for = 

od iy irrgation of crops. The lowest, 42 gallons, was in the evirons of the village 

7 j _ of Patchogue, a relatively undeveloped community. The average | use per 

capita for all the water-supply plants tabulated was 76 gallons. 
= | —_ More than three- -eights of the nite of Suffolk County depend entirely upon 


| mestic wells probably will be wer 


— FF increase the water use per capita. Consequently, estimates of the future 


water-supply needs of Suffolk based on present-day use capita may 
be seriously it in error. | 


4 
= Withdrawals from private industrial wells supplied about one-third ofthe | 
nto total ground water pumped in Suffolk County in 1955. These withdrawals 
: ; largely took place in the western half of Suffolk County, being used by hos- al 4 


~pitals, aircraft factories, food-processing plants, dairies, laundries, light 


manufacturing plants, storés, offices, and public utilities other than ‘water 


cooling including 2 air conditioning. As cooling does not alter the chemic cal 
content of the water, the State conservation law requires that the spent water 
be recharged to the source formation. . This is commonly done through wells, ; 


diffusion wolle, and excavated pits, , called recharge basins. 


About two -thirds of the water pumped for industrial purposes i is used for ix ; 
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Greenport Village of Greenport | 166.1 


day period starting in the middle of June. . 


Since the end of World War II, acreage irrigated by water in Suffolk 
‘ County | has ‘more than quadrupled. . The most intensively irrigated area of a 
ie Suffolk County is the northeastern corner, north of the town of Riverhead and 
‘ Pigure 4 shows a typical irrigation-well installation. The pump which is - 


| mounted over the well on a concrete base, is driven Gisect 5 
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7 Ss _ The third important use of ground water in Suffolk County is for irrigation. tio 
1955, one-quarter of the total ground-water withdrawals was used to irri- 
4 a _ gate about 56 square miles of intensively farmed land, about one-third of all : re 


- 


Fig. 4, Pump | used ‘Supply irrigation water 
‘mows the by well is —— is pumped through 
_ portable aluminum pipe to oscillating sprinklers which spray it in a circular — 
pattern, Sprinklers differ widely in capacity and design. ‘The type shown in 
= the figure distributes about 15 gallons per minute. Another oil can dis- i. 


400 gallons per | minute through o one Sprinkler | head, 


Some farmers irrigate when the soil no longer cakes upon being eed. 
Others follow the recommendations of the Long Island Research Farm at a 
Baiting Hollow and supplement the rainfall with enough water to insure the 

a application of 1 inch of water per week to the land. Still others, who have few 
- pieces of equipment and are short of manpower, adjust their wie pro- 
_ No data are available but it is thought that a large percentage of the irriga- 
5 tion water applied to crops in Suffolk County is evaporated from soil and al 
5 surfaces, At times excessive amounts of water may be applied to the . 
‘resulting in the return of part of the water to the water table. In general, a 
‘a however, probably only a small part of the water pumped for irrigation gg 
returned to the water-bearing deposits. 
_ The rapid —oreaee in irrigation in the last 5 years is evidence of the re- 
mite obtained. The farmer not only protects his crop investment during ary 
years, but also” increases the yield per acre in years of normal rainfall. 
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. “Suffolk County are sea-water contamination, and corrosion due to the use of 
4 untreated low pH water, Once full- ‘scale development of the ground-water 


eo reservoir begins, additional problems m may become important. 


Sea-water contamination has been observed only in isolated in near - 
shore ‘areas or in other areas underlain by salty water at shallow depth 
s : (Hoffman and Spiegel, , 1958). ‘These a areas include p parts of the North and 
_ South Forks, off-shore barrier beaches such as Westhampton Beach and Fire 
‘Island, and a | narrow zone, of undetermined \ width, adjacent to and i paralleling 


the: shoreline of Suffolk County, where the water table is 3 only a a few feet mae 


™ 

wnt 
— cha 

— res 
— 

i= 
wid 
2 
4 Ast 
pret 

ai: 

— 

war 

—  slig 

q 

__eritical contamination in on 3,000, near the eastern tip asl 


¥ the North Fork, utilizes about 130 to 160 ‘million gallons of of water a 
_ At 3 of the 4 well installations, the water pumped has chloride voeumaaes 

markedly above normal. During times of peak demand the water pumped at 2 ol 

‘a these stations has chloride concentrations of more than 400 parts per mil- 

- ton, considerably above the maximum of 250 parts per million recommended 4 
_ by the U.S. Public Health Service. _ Only be adept combination of this water 

with water of low chloride concentration from the fourth pumping station can 
a potable water supply be delivered to the user, 

a previously shown, much of the ground water of Suffolk County is cor-— 

. rosive. Large water-supply installations can cope with this problem eco- 


_-Romically, however, corrosion treatment for individual well systems is | not 


A 
* ad of these systems must be replaced ultimately. _ Although no specific dollar | 
3 fe value can be an the resultant financial- loss can be considerable. em be 
arising from increased ground-water utilization in 
| folk County can only be conjectured. Some of the more serious of these are: 
(1) more extensive sea- -water contamination of the ground-water reservoir; 


(2) contamination of the ground water by industrial wastes, cesspool effluent, 
and fertilizers; (3) increases in ground-water temperature due to on 

recharge; and (4) permanent depletion of ground- -water 
‘Suffolk County should be extensively sewered, sia Mad 
_ Heavy ground-water withdrawals, particularly in shoreline areas, can 


cause extensive sea-water contamination of the ground-water reservoir. 


is s underlain by even heavy development may 
— cause encroachment, On the other hand, extremely heavy withdrawals in the | 
= inland portions of Suffolk County may produce no measurable en- 4 

Ground- water contamination can also result from the recharge of at 
| In future years, , without remedial measures, wastes of all kinds may be re- 
charged to the water-bearing sands faster than they can be discharged. | AS a 
result the concentrations of certain compounds will increase. Persistent _ 


ae from this source ultimately might create local, if not sean ¥* 


= 


— 


i. = the present time, contamination by high-chloride fertilizers is indicated is 
_ in the intensely cultivated areas of Suffolk County (Hoffman and Spiegel, 1958). 
“ _ As the history of this contamination is largely unknown, very little canbe 


predicted concerning its progress in future years. _ However, the | concentra- oe 
tA tion of nitrates and chlorides in the. ground water possibly 1 may y increase —_— 
ficiently to render the water unpalatable or to cause organic disorder. 

ead Changes in the temperature of ground water can be as 3 important as chemi- . 


cal contamination. recharge of spent cooling and air- -conditioning water 

g _ _ Warms up the ground water in the vicinity. The recharge of small volumes ot 

4 K slightly warmed water creates no problem, whereas the recharge of large onal 

= volumes of high- -temperature wa water does, and in intensively developed indus- 7 
he 


~ trial areas this condition could result in costly cooling problems as happened © 
at one locality in western ‘Long Island, N. where well- -water 
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and of ~water storage in Suffolk County 
may r result from extensive sewering. _ ‘The only sewers in Suffolk County at She. 
the present time are the small systems of the near- -shore villages of ir ge 
Huntington, Northport, Ocean Beach, Patchogue, Port Jefferson, and 
Riverhead. So, probably 80 percent of the v water pumped from public, ; Private 
_ and industrial wells is returned to the ground, If the area served by evwers 
expands, and if storm water is discharged directly into the ocean, replenish-— 

_ ment previously obtained from sanitary wastes and storm water will be con- 6 
4 On the other hand, withdrawals from the ground-water 
& reservoir will not lessen and the increased net draft on the ground-water will 
a eat cause a decrease in ground-water storage and a decline of the water table. Care 

One consequence of appreciable -water- -level decline would be the need to 

= some of the existing shallow wells for further use. Another would be | 
an increase in pumping costs incurred by the additional lift required to raise 

- water to land surface. Perhaps of greatest consequence is the possibility 
_ Of a landward movement of sea water that would cause the complete — 

_ An example of the problems that can arise from intensive development is is ; 
in ‘Kings County, the westernmost county of Long Island, New York. Here, ait 

_ the population density is more than 36,000 people per square mile. - Ground 


water recharge is greatly reduced by paved surfaces and the sanitary and 
storm sewers that drain the area. Heavy pumping by the New York Water 
Service Corporation, near tt the geographical center of the county, lowered ; 
water levels at some places to more than 20 feet below sea level. The result- 
ing sea-water encroachment raised chloride concentrations in the adjacent 


oote water to well above 1, 000 — per million, and the system was shut 4 


problems described d above cé can be costly. ‘Possible remedial meas-— 
ures might include: e: (1) good water budgeting; (2) the proper location of — aa! 
ground-water withdrawals; (3) the removal of wastes by sewers; (4) constant as 
supervision and evaluation of | ground -water conditions; and (5) foresight in nf: 
a Good water budgeting serves to increase ground-water replenishment — vi 
while decreasing ground-water losses, If If present growth trends persist, at 
some time in Suffolk County’ s future, seaward disposal of sewered wastes 
4 will probably be required. Consequently, concurrent with sewer 11 installation, © 
“Measures to offset ground -water depletion should be considered. These 
a might include properly located storm-water recharge basins and, if ac- HS e 
ceptable, sanitary-waste leaching basins. _ Recharge | near the shore of Ves 
- selected part of the sewage would help to buffer sea-water encroachment. — 
_ Another budgeting measure to offset ground-water depletion would be to a 
_ duce the natural ground-water discharge by streamflow, underflow, and ay 
ged Wells near the shore or next to pumped at appropri- 


Careful selection of well-sites and regulation of 


drawals would help to minimize sea-water tides one hand, well 


the | same time make the land available for development. 
aa 


+ 


3 
| 
3 sites, 
Some 
socio 
presi 
| cals 
: 
unde: 
ultim 
| 
authe 
— Cour 
— Ville 
— 
— 
— 
— 
4 a forr 
and 
— 
Vg 
— 
— Yor 
“4 
a i 
— 
= 
a 


& 4 


- 
7 | 


GROUND 


pp’ should be located a as far from the shoreline as possible. On the other 
hand, to intercept natural discharge, they should be placed near the shore. — 
Some adjustment is possible in the type of installation constructed. 
complex interrelationships of the geologic, hydrologic, hydraulic, and 
socio -economic factors involved in the development of Suffolk County’s a. 
- ground-water reservoir requires careful evaluation of the problem. — At the 
present time part of this evaluation is being undertaken by the U. S. Geologi-_ 
cal Survey. In with local agencies, the Geological Survey is ap- 
_ praising the occurrence, , movement, quality, and quantity of the ground water 


ies Suffolk County, in addition to investigating conditions which might | 
ultimately threaten its quantity or quality. 
_ The New York State Water Power and Control Commission, which has om { 
authority to regulate large-scale uses, except agricultural, supervises reser- 

Ve voir development. This regulatory body licenses well drillers, and deter- a 
mines policies concerning ground- -water and artificial 
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4 _ DESIGN METHODS FOR FLOW IN ROUGH CONDUITS* 


Henry M. Morris,* F. AS ASCE > 


ey 


Design curves and methods are for determining fri friction factors 
ft for turbulent flov flow in closed conduits and tranquil open channel flow. 1 The ~ 
- methods are » physically realistic, being based on five distinct regimes col 


o of the most basic and frequently occurring engineering is 
that of the conveyance of water or other fluid by means of some sort of flow 
channel. The channel | may | be either closed or open, the basic mechanics 
being | essentially the same. Usually the conduit lining is rough 
ac _ Despite the commonplace nature of ‘this problem, , design m methods for deal- ete 
ing with it are still largely e: empirical. Formulas such as those of Manning, eS ‘a 
Kutter, Scobey, Hazen- Williams and many others are illustrations of this ae 
fact. Each of these formulas is dimensionally incorrect and has a relatively | e. 4 : 
limited range of physical application. Each applies only to the turbulent flow 
of water, at relatively high Reynolds Numbers, and assumes uniform flow con-_ = 
ditions, But their most serious defect is s the wholly empirical nature of the 
Provision made for the effect of boundary roughness. It has often been shown | 
these so-called “roughness coefficients” not only vary tremendously 
3 “the specific nature of the conduit surface, but frequently vary also with the 
ynolds Number, hydraulic radius, and omer factors. 
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 ‘Itis well- -known that dimensional requirements establish | correct f 


function to be of the form known as the Dare e 
y 
where H,/L represents the unit expenditure of overcome 


diameter o of conduit. The in the us use of this ‘simple and universal 7 
: formula is of course in the evaluation of the dimensionless friction factor ‘ 
Which can be shown to be, in (assuming uniform flow), | given 


The di dimensionless force 1 ratios, Nr, NF, Nw) NM represent respectively the 
‘2 aa influence o on n the flow structure of the forces of viscosity, gravity, surface — ae 
= and fluid compressibility. For most problems in civil engineering _ 


—— is the ratio of “se head to wseeelet 


hydraulics, the last two have negligible effect and can be eliminated from the 
function. Gravity has no direct effect in closed conduit flow and is often 


negligible in problems of tranquil open-channel flow. It will be assumed 
herein that this is the case, so that the following discussion will not apply ba 


; “a quantitatively to rapid flow in open channels. The Reynolds Number, Nr, 4 
however, must nearly always be considered, since every flow phenomenon > 
iy 7 necessarily involves an expenditure of flow energy in friction. The other © 


terms in the function reflect the influence of the wall roughness elements on 
the flow and turbulence patterns. Thus h, A, s, r, represent the radial height, 
eS longitudinal spacing, peripheral spacing a and radius of curvature of the rough- 
‘ness elements. As many similar additional terms should be included in the wd 


function as necessary to describe completely th the geometry of the roughness Bing 

S Because of the complex nature of the friction factor function, various at- mad 


7 tempts have been made to group all the roughness dimensions together into a 
singe representative dimension for simplicity. . Since, however, each of the — 
various dimensions may have a significant effect on flow geometry, this rep- 
a resentative dimension usually cannot be any actual dimension, and must there- 
fore be determined empirically. _ It thus becomes little more than another 
roughness coefficient. Furthermore there is no real assurance 


; ~ :meter” as the measure of roughness is a case in point. concept, 
based on the well-known Nikuradse tests(1) on sand- coated pipes, has been — 
_ highly popularized by the inclusion of the “Moody curves” in nearly ili —n 


a ie These | curves, which give the friction factor as a function of the Reynolds — 


_ Number and the ratio of equivalent sand roughness to pipe diameter, oil 


— are implicit in this procedure, as follows: 


r; 1) The Moody | curves and the empirical values of equivalent sand roughness 
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am turn is merely ‘that of an artificial curve asymptotic to the Von meee 
‘ a Nikuradse so-called “smooth” and “rough” pipe curves. This is a Lie 
= transition function, showing the friction factor to be a decreasing func- @ 
- tion of the Reynolds Number. The artificiality of this device is evident, 4 
: since the actual experimental transition function of Nikuradse, on ‘whose 
results this procedure is based, is exactly opposite in character to the come 
- ‘Colebrook function, showing the friction factor as primarily an increas- 
_ ing function of Reynolds Number. _If there were really anything “ Fequiv- 
alent” about Nikuradse’s sand roughness, it should at least show the 
_ same form of functional relationship as obtained on the commercial 
2) The Colebrook-Moody curves make no provision for conduit surfaces 
_ producing a rising friction factor- -Reynolds Number relation (such as 
- corrugated metal pipes(4) or flumes, sand-lined channels and conduits, = 
) or those with a fully horizontal characteristic(5) (such as conduit 


He 


3) Even for surfaces which yield descending curves, there are e many 
indicating that the Colebrook equation is not an 1 adequate descrip-_ 
tion of this relation. Some tests on large conduits show that friction 
nt a factors tend to decrease up to much higher values of Reynolds Number “J 
than predicted by the Colebrook relation.(7) Certain tests have yielded 
x ‘systematic variations in empirical values of equivalent sand roughness 4 
a - for a given surface, (8) proving it to be unreliable as a means of pre- | a 
iz x dicting friction factor unless these variations are taken into account he y 


= ‘Various data are available(9) refuting the Colebrook Equation implica- — 
r * tion that friction factors for a given surface type always decrease with 


i- 


Ss Perhaps the most serious objection to the use of the “equivalent — 


ness” ’ concept is that its superficial rationality tends to give the de- ~ 

ay signer an unwarranted confidence in the accuracy of its predictions. = 

The large number of rougliness dimensions influencing the flow turbu- 


1 


7 lence indicates the impossibility of representing the element by any — wai ° 
dimension, - especially some non-existent “equivalent” dimen- 
sion, - and the several types of possible f-Np relationships indicate 
a the impossibility of using any single function to define this relationship. 3 + 


ve turbulence will show that there are several distinct regimes of turbulent flow a 
possible in conduits, each with its own peculiar characteristics. Turbulent 

- flow is characterized by the continual generation and dissipation of vortices, = 
with most of the energy loss occurring during the process of generation, the 

: remainder in the process of viscous dissipation. In general, the generation of 3 


vorticity is a viscous phenomenon resulting from fluid moving at relatively P ae 
high velocities in contact with fluid moving at relatively low velocities. The _ 
x phenomenon typically occurs either at a laminar boundary layer or at a zone — Si 


; bh either case, the vorticity so generated washes out into the main 1 body of ‘4 
f=. flow. The processes of mass and momentum | transfer, as well as viscous aa 
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r attrition, combine to reduce the turbulence structure toa certain typical pat - 
tern in the interior region of the conduit, regardless of the method by which © 


he original vorticity was generated. 
This *‘normal” turbulence can be statistically, yielding tl the well- 


10wn velocity distribution equation: 


“In this equat ion, vy is the waineiiee at Manone: from the he wall, Ve i is the cen- 


ter ‘velocity, at distance To foom the wall and Vs is s the | “shear velocity,” /de~- 


turbulence constant characteristic of tie “normal turbulence.” Equation 3), 
Fi a if plotted on semi i -logarithmic paper, , yields a straight line, of ‘slope 1/k. I a 
is assumed, of course, that the section under discussion is sufficiently re- 
moved from the inlet or other disturbance to permit full development of the we 


_ However, normal turbulence will usually prevail only in the central 1 re- - , 
gions of the conduit. Near the wall, the velocity distribution will be material- 
a 74 ly affected by the nature of the vortex -generating mechanism. If the wall is 
as smooth, with a laminar boundary layer, the very small-scale vorticity being 
d| J generated, together with the transitional effects from laminar to turbulent 
flow, causes a steepening of the » velocity gradient for some distance from the 
‘On the other hand, if the wall has roughness elements the 


laminar film, large-scale vortices are continually generated in the separation 

= zones behind each element. ‘This is manifested by a relative flattening of the 

velocity gradient for a region near the wall. These various velocity distribu - 
phenomena a1 are illustrated i schematically in Figure 1. 


effects would be combined, It becomes obvious, then, that the longitudinal _ 
_ frequency of vortex- -generating roughness elements is of determinative sig- _ 
nificance with reference to the regime of turbulence that will prevail. With aed 
these preliminary considerations in mind, several — of turbulent flow _ 


Gy cage 
in conduits may be be considered. 


If the conduit wall is smooth, so that the turbulence - 
is solely | one of minute vorticity shed from the laminar boundary layer, then — 
a the velocity distribution and friction factor are independent of all wall rough - : 
a ness dimensions. They depend rather on the thickness of the boundary layer, © 


which is a decreasing function of Reynolds Number. The 


The constant A is empirically determined to be 5.5 the turbulence 
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METHODS 


BOUNDARY 
10 12 14 16 


riv or d its. The 
“equation fo for open channels( 10) i ist 


However, equation (5) m may be used for open channels without serious error, 


turbulence is defined as that that characterized by t the universal 
_ turbulence constant k. SS Normal turbulent: flow would therefore involve ay 


— 

| 

ie 


=) 


could be no wall zone e of whee ‘mb- -normal or abnormal turbulence. Vortices. 
 gueevenad at the wall are almost immediately broken up into patterns statis- 
_ tically equivalent to the normal turbulence of the central region. ae pia 


7 The velocity distribution equation for normal turbulent flow may be written: 


in the case sand as used by In 
mal turbulent flow the other roughness dimensions are irrelevant, except that 
_ the crests of the elements determine the physical limits within which the 


= distribution equation is applicable. _ Whatever shape the elements may 


_ Experimental work has yielded a value of about 8. 5 for the constant B, for 
va rious types of roughness elements. " The value ‘of about 7.8 seems somewhat 
tter for two- dimensional closed flow or rectangular open 


rresponding factor equation(11) is; 


. has no effect on the turbulence pattern, which is fixed in terms of the — 


Here 


fined in terms of the actual roughness not some equivalent 


an ordinarily smooth conduit surface is with occasional 
- isolated roughness elements, the o over - -all friction factor will be that due to. 
= friction drag at the laminar mee layer plus that due to the form Greg 


= 
"where fs is the smooth turbulence friction factor at the 
Number, E an element ch characteristic defined as 
the pex yer ipheral of ther: 
roughness | elements to total conduit | wetted perimeter). ro is the re relative — 
: roughness height, ro/x has a already been termed the relative roughness yf xii 
Spacing, and is ‘the roughness index. The coefficient Cp is the drag 


is the longitudinal spacing of the roughness elements, which is the 
roughness dimension characterizing the frequency of wall vorticity sources 
4 
a 
e case ofa non- 
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“a value of E can be computed for each repeating type of roughness aes - ; 
"ment on the surface. . The individual values are > added to ti the sum effect. z 


similar to that for smooth flow, with the friction factor computed as the ot Pe 
smooth turbulence friction factor multiplied by a number greater than unity, i 
the value of the number being (1. + 67. 2 LYE). The friction factor for semi- | 
e smooth turbulent flow therefore normally decreases with increasing Reynolds ~ 
the wall roughness are e sufficiently close together, the wake be- 

hind each may extend to or nearly to the next element. _ There is then essen- ie 

_ tially no part of the wall over which a laminar boundary layer exists. © Fur- 3A 
B thermore the vortex generation and dissipation phenomena associated with — + 
each wake will interfere with those at the adjacent elements, so that the ne in= 


dividual effects are not additive as in the case of semi- smooth flow. cine ie 


{35 
EMI= SMOOTH TURBULENT] ae 
| 
| 
4 
| 
_ Wall of abnormally intense turbulence and mixing ribution 


as 
will he. normal in the central regions, but the a average slope near will 


n for this type c eof 2s can be shown(13) to be 


‘ ing in this zone. This type of flow may be called Lt -turbulent flow, and i is ri 


‘The friction factor equatio 


rs 


Tae 


“The function 7) can be shown to decrease with increasing Reynolds Numbe 


ren 
NR is sa wall all Reynolds Number,” to 


‘The fu function may be regar asa transition function whose magnitude 
4 


to increase 2 the increasing Reynolds Number, whereas for both smooth and arty 

semi-smooth flow the friction factor decreases with Reynolds Number. 

_ seems likely that for surfaces which will yield hyper - turbulent, and ultimatel ay 

normal flow, there is a brief regime of smooth and then semi-smooth er a. 3 

before hyper-turbulent flow begins. This sequence is evidenced by the _ 

familiar dip - -and-rise of the Nikuradse rough pipe transition function. 

‘Hyper -turbulent flow is produced over such surfaces as as corrugated metal, 

Bey sand- -coatings and other surfaces of considerable roughness. The exact form 

transition n function depends on t the > shape of the wake zones behind the 
elements and therefore on the form of roughness elements, and must thus far 
be determined experimentally. It is very significant, however, that for suf- _ 4 

a ; 7 a high Reynolds Numbers, each transition | function approaches the 7 
normal turbulent flow equation, for which the friction factor depends coed ie 

on the relative roughness spacing. 

a still another type of flow must occur over surfaces ‘composed either of i 
small depressions or of roughness elements spaced socloselyastoforma 
more or less smooth pseudo-wall composed of the element crests and the en- “i 
of dead fluid. Within these pockets or depressions will + 


_ tices will be maintained through transmission of ne stress from the flow- ~£ 
: bs ing fluid at their upper limbs. In addition, small-scale vorticity will be gen 4 
erated continuously along the -wall somewhat analagously to the 
7 process in smooth turbulent flow. ‘Fig. illustrates the = 
energy expenditure is partially to 
; . ticity at the pseudo-wall and partly to maintain the stable depression vor- 


tices. An + oad anerame for the friction factor in this | type of flow 
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Ih (12), | jis ‘the longitudinal dimension of the is 
a velocity near the pseudo- -wall, and Se. is the coefficient less than unity such 
that CwVy is the peripheral velocity of the depression vortex. Other terms 
4 have been defined previously. Examination of the quasi-smooth flow Rey 


reveals that the friction factor will normally | decrease with incnanning, Rey- 


£44004 


Design | Curves. for Flow 


In order to to facilitate design calculations as as 
have been plotted for each of the various turbulent flow regimes. = These are ; 
= the same general type as the familiar Moody curves, giving friction factor . 
_ as a function of Reynolds Number, and are used in exactly the same way. The 
™ only difference is in the roughness parameters. As discussed previously, — oa 
is incorrect to attempt to. lump all roughness effects together in terms of any a 
z kind of equivalent dimension, Each flow regime requires both its own rough- a 
ness parameter and friction characteristic curve. 
= Fig. 5 contains the curves for semi-smooth turbulent flow, with the lower- 
_ most curve being that for smooth flow. The curves are plots of equation @), 
for parametric values of the roughness function J E defined in equation (10) 
- Values of ‘drag coefficient. Cp (ordinarily assumed constant with — 
_ accuracy for design purposes, at usual values of Np) to use in evaluating “— ; 
for each repeating roughness" element are tabulated in Figure | 5-a. Use 
these curves of course requires s that the element height be great 


vik _ The thickness of this boundary layer i is a decreasing f function* * of yy 
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SEMI-§ “SMOOTH FLOW EQUATION: 


FLOW FRICTION FACTOR RAT 

“PERIPHERAL ROUGHNESS __ RATIO 


"CONDUIT WETTED PERIMETER 


"AGING a 


ELLIPTICAL 


a = FOR USE WITH SEMI-SMOOTH FLOW 


a 
| 
STRIP ROUGHNESS 
— 
q 
— 
aa 


Reynolds Number. This function is plotted on Figure 5- 
roughness height is significantly smaller than 6, itcanbe assumedthat 
smooth flow will prevail. For simplicity, and for design purposes, no - 3 
Bi. has been made to delineate i a transition function from smooth to semi- 
smooth flow. . For particular elements, this could in theory be done by allow-_ 
ing for varying values of h and Cp at low Reynolds Numbers, until the ; 
points: on the crests of the elements become stabilized. 
_ The quasi- -smooth function is shown in Fig. 6, again with the lower- most 
“curve being that for smooth flow. These curves are Plots of equation (12), 
with the roughness parameter as defined therein. ratio Vw/V may 
range from 3/8 to 6/8, perhaps averaging 2/3. The coefficient Cw is yet =- 
known but would perhaps range from about 1/2 to nearly unity. > | . 
a} — (Cyvy/V)3 would therefore lie within the range of, say, 0.01 to 0.40. In the 4 
ie _ absence of better information, a value of 0.10 would be probably reasonable : 
7 and conservative. The value of 0.05 checks well with such observational data 
as are currently available for this type of flow. 
Since the friction factor in hyper -turbulent flow depends on both rela- 
_ tive roughness spacing rg/ and the form of elements, a different set of 
curves: would be required for each form of element. On ‘Figures 7, 8, and 9 = 
“curves have been plotted for the most common on types of rough surfaces. 
4 ducing this type of flow, as well as its limiting state of fully normal turbulent 
; flow. . These surfaces are those containing regular | corrugations, sharp- edged at 2 
circumferential and densely -packed spot roughnesses, such sand 


om Each set of curves is s plotted with the relative roughness spacing to/h as 


: parameter. > ‘If the effect of some other form of roughness is needed, a i ae 


_ judicious extrapolation or ces ag of the results obtained from ‘these 


curves: (which have been shown\ 15) to correlate experimental results obtained 


from n many different sources) 1 might be > employed. it ‘may be be noted too that Pol 
the curves, since they are experimental curves, define the entire transition | 
_ from ‘smooth flow to normal turbulent flow. As noted before, this includes: a 


= of semi-smooth flow as the hyper- -turbulent. regime 


. othe e curves of Figures 5 through 9 will permit the on solution of any ordinary 
; _ problem in conduit flow (in which elastic, surface tension and gravitational 
effects can be neglected), provided only that the roughness geometry and ell 
— are known. However, the problem of discriminating between semi- —, ~ 
a smooth, quasi- -smooth and hyper- -turbulent flow has not yet been discussed, = 
- Quasi- -smooth flow will occur over roughness depressions in which stable — 
- vortices « can be maintained, This ager that these depressions must be at © 
least as deep a as s they are long. Otherwise they will develop either 
N : definitive, that at sufficiently high Reynolds Numbers normal turbulent flow - 
may under some conditions change to quasi-smooth flow. This i is indicated z A 
_ when a rising or horizontal f-Np curve begins to drop. 2. Until more informa- 
- tion becomes available on this condition, it is conservative to ignore it _ ae ¥ 
ee The problem of determining whether, for a given yen surface, semi-smooth or | ro 


_ hyper-turbulent flow will occur is more difficult. For a oa type of rough- ees 


ness element, the spacing of cloments is bane critical factor. 
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that spacing for which the friction factor attains its maximum value. -. Joe 
3 spacings exceeding this critical value, | causing semi- -smooth flow, “equation (9) 
: shows friction factor to decrease as the spacing increases. For closer spac- - ar 
_ ing, producing hyper-turbulent or normal flow, equations (11) and (8) show t the oe 
== factor to decrease as the spacing decreases. = — 
The critical ‘Spacing (or perhaps another factor) for a given type o of ele- 4 
_ ment could therefore be determined by equating the friction 1 factor expres- aos? 


sions for semi- -smooth and hyper -turbulent Thus: 


| 


; between the two regimes, for : a given Reynolds Number and conduit would be 


, 
Solution of equation (13) for the spacing A, or whatever factor is unknown, _ 

determines the boundary between semi-smooth and hyper -turbulent flow for 4 

the given conditions. The presence of the function @ makes such solution 


quite difficult. Asa reasonable approximation, assume fg 0.01 and dd=0. 
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Applications in Practical Hydraulic Analysis" 


‘The application of these concepts and  denign curves to practical flow prob- 


lems requires knowledge of the geometry of the surface cs omanrane 


produce hyper-turbulent flow, and at sufficiently high Reynolds Numbers, | - 

i: Materials commonly used for water and ‘sewer } pipe, and for gas and oil te 
transmission - including concrete, cast iron, clay, welded steel, etc., - are, : 


least when new, essentially smooth, except for isolated roughness 


ments and for the joints. If rivets, bolts, etc., are present, their effects are 
easily computed by the semi- -smooth flow equation. The joints may cause 
either a protuberance ora depression, and the effects may be calculated by 
the semi-smooth or quasi-smooth relations accordingly. “Although the exact 
wa geometry of the elements often may not be known, it is believed that the de- _ 
a Signer” could estimate these factors with at least as much confidence and ac- si 
m8 ‘curacy asa a roughness coefficient or equivalent sand roughness. . Materials 
is of this type will usually produce either semi- -smooth or quasi-smooth flow, _ 
the appropriate curves can be used accordingly. 
‘Surfaces of variable or irregular roughness are of course more difficult — 
to evaluate. However, the principles of the various turbulence regimes must — 
_ still be valid, ‘so that : friction factors for * such surfaces. could still be esti- 
mated by judicious application and extension of these principles. 
If, for example the surface consists of a combination of several types of 
isolated protuberances and depressions, the ‘semi-smooth and 
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equations could be as follows: 


_ _ If the elements are close enough together to produce hyper-turbulent c or = 


- normal flow, the hyper-turbulence equation must be modified to allow for _ 


- average value may be used with sufficient accuracy, since it is basically the 

_ number of turbulence sources that determines the wall zone turbulence pat- 

_ terns. If both the spacing and roughness form are variable, the average bres 
value of spacing may again be used in the hyper - -turbulent flow equation, but — 


— of various forms and spacings. If only the spacing is variable, its. = 


™ When occasional elements project much farther into the flow than oe 
s effect is that of superimposing isolated wakes and drag forces on the typical } 
 hyper- turbulent flow near the smaller elements. It is reasonable in this im 
‘ situation to compute the friction factor due to the hyper -turbulent flow near a 
the ordinary elements, and then add increments of friction factor corres Be 
?P ponding to the elements. ‘The resulting would 
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equation, using the of larger elements on 
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the velocity | near the wall which impinges on the projecting 
74 about 0.7 of the a average velocity in the cross-section. The element height h- 


_ must be measured from the crests of the smaller adjacent elements. 
In the case of a conduit with part of its perimeter very rough and part See 
porns (such a asa a corrugated p pipe with paved invert, or a pipe with longitudi - . 
nal rows of closely-spaced rivet projections) it is obvious that the over-all 


a factor will be intermediate between that for a smooth flow and that - 


for a hyper-turbulent flow corresponding to the rough section of the surface. — = 
_ Most of the energy expenditure due to friction in a flow occurs in the process 


of vortex generation at the conduit surface. ‘Thus it is reasonable to take the - 


proportion of rough surface area to smooth surface area as the basis for d 


_ termining the bulk friction factor. The resulting equation is: 


2lo — + 1.75 + 

“where P, and p are the smooth and rough segments of sealtienenes: > 
In the relatively rare case of a surface composed of “random” ” roughness — 

_ elements, of various shapes, sizes, and spacings, the tendency is for the > wild 5 


friction factor Reynolds Number curve to be horizontal throughout the en 
ie tire turbulent range. That is, such a statistically random roughness distribu- 


tion produces a statistically random vorticity pattern near the wall, and S 
therefore normal turbulent flow, even at low Reynolds Numbers. The friction 


factor for such a surface can probably be estimated by the normal turbulent os 
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- 3 The concepts, equations, and curves presented in this paper provide a de- | 


<a sign and analysis tool for flow problems which is substantially r retional in na- 
ture, as compared with the largely empirical methods in common use. It is : “ 
somewhat more difficult to use than previous methods, a defect inherent in | & 
_ any attempt to place the complexities of hydraulic phenomena on a more ra- . os 
" tional design basis. However, this difficulty is ; materially alleviated through | re 
4 use of the design curves included herein. 
The slightly increased effort required of the designer is, , moreover, more 5 
{ oe. justified if his work can thereby be placed on a sounder scientific footing 
4 ry than is possible near present procedures, which necessarily involve an im- a 
portant element of uncertainty in the choice o of a “roughness coefficient,” ; 
; “equivalent sand diameter” or other roughness measure, as well 2 as ‘uncer- 
“tainty in the variability of this measure. | 
‘The advantages and potentialities of the concepts of this paper, for solution _ 
of practical problems in both open and | closed conduit flow, would appear to = . 
_ justify intensive experimental study for verification purposes, It has already a 
; _ been shown(16) that the concepts correlate data obtained from many different us 
e sources, On many types of roughnesses, in both open and pressure channels. 
my It is believed that this is also substantially true for ‘more recent experimental 


g refining a> senna te and design methods based on them is s very desirable. _ x 
Ki “Wa Until such further studies become available, it is believed that judicious 9 
use of the equations and curves presented | herein will suffice even at present | 
3 to give ve design results at least as reliable, mn often more than any 
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=a constant in the smooth flow velocity distribution equation, | taken as 


*~ = a constant in the hyper - -turbulent flow y velocity distribution equation, 
ees eee as 8.5 for axial and 7.8 for two-dimensional flows. hiss a 


= discharge coefficient for surface roughness Saas . 


= coefficient less than unity such that Cy Vy is the peripheral 


roughness clement function = Cp DP 
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eight of roughness elemen sepa 


‘= von Karman turbulence constant, taken as 
= length of reach in 
Number (used as Dv or 


r to = ‘Mach Number 


lready 
P = perimeter of conduit 


- total peripheral length of roughness elements at a secti 


=e uivalent radius of conduit 
radius of rounding of roughness element 


= peripheral of ‘smooth wall surface 


_= shear velocity, whose square is the — of wall da ian -stress 
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= (turbulent) close to the wall 


depression. vortex function in 
additive element in n hyper 


-center-line longitudinal spacing of 
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] Numerous experimental curves of these types are shown in yee i 


Monograph 7, “Friction Factors for Large Conduits Flowing Full, 
N. Bradley and L. R. R, Thompson, S. Bureau of Reclamation, 
“High Velocity Tests ina _Penstock,” by Maxwell F. Burke, Transactions 
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- ae: Pipes,” by L. G. Straub and H. M. Morris, Technical Paper No. 3, Series | 
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“Hydraulic ‘Tests on Concrete Culvert Pipes,” - LL. G. Straub and H. M. 

‘34 Morris, Technical Paper No. 4, Series B, St. Anthony ‘Falls lls Hydraulic 

Univ. of Minnesota, Minneapolis, Minn., 1950. of 
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American Geophysical Union, , Vol. 30, 1949. The equation given is 
Powells Equation (7), with the Chezy C replaced by 8g/f 
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_ GROUND WATER FLOW IN THE NETHERLANDS COASTAL DUNES 


SYNOPSIS 

mathematical analysis of the ground water in that. portion 
“the Netherlands coastal dunes which furnish a water supply to Amsterdam is her 
"presented. > The problem is treated two-dimensionally with a fresh-water lens 

_ above and a fresh-water pocket below a clay layer under the dunes. Salt water — 
underlies the entire area. _ Beginning with the basic equations describing the — 
ground water situation, assumptions and boundary conditions necessary to ob- Pa, 
solutions are described. Results give water movements and 


— 


isa portion of the dunes along the North Sea coast 

_ Southward from Zandvoort for a length of 10 kilometers. Beginning in 1853 

water was collected from a series of shallow drainage ditches; subsequent ¥ 

expanded the ditch system and added deep wells. 
Fig. 1 shows a cross- -section of the area perpendicular to the coas coast. The. 3 

> dunes parallel the coast, are about 4 kilometers wide, and are mostly 5 to 10 I 
meters above mean sea level with maximum elevations of 30 meters. ‘The ; 
average height of the ground water table is only two meters above mean sea . 
level with the highest level at about six meters. | ‘The dunes are bounded on the | es 
west by the North Sea and on the east by a 4-kilometer wide strip of flat land .< 

ig with a ground water table varying from +1.0 to -0.6 meter mean sea level. 4 
East of this strip lies the Haarlem Lake Polder with a ground water level a ii r 
meters below mean sea level. 

es Note: Discussion open until December 1, 1959. To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE, Paper 2082 is ‘oe 


a" part of the copyrighted Journal of the Hydraulics Division, Proceedings of the - a . em! 
American Society of Civil Engineers, Vol. i. Me. HY ‘July, 1968. ag 
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DISTANCE FROM THE COAST, KILOMETERS 
ha fit 
SECTION OF DUNES PERPENDICULAR TO THE COAST SHOWING 
WATER DISTRIBUTION IN. 1850 ( AFTER 

— large number of wells and test holes have enabled the geologic structure 
of the dunes to be defined. In terms of ground water movement four strata can 


a be distinguished. The dunes are an aeolian formation and from ‘ground surface 


@ 
= to 13 meters below mean sea level consist of a uniform fine Holocene a 


A much less pervious stratum and peat occupies a | layer 


‘entends to 160 meters below mean sea 2 level. _ The final stratum, extending 
downward from 160 meters below mean sea level, is an Icenian (Pliocene) 
most of the country - is s underlain by salt water. Precipitation infiltrates into 
the dunes to form a lens of fresh-water above the clay layer, while a portion | 
of this water percolates through the clay layer to form a fresh-water pocket 3 
floating on salt water. The distribution of fresh and salt waters in 1850, be- 
fore withdrawals for water supply were initiated, is indicated in Fig. 1. __ 
A considerable body of literature relating to sea water intrusion and fresh- 
_ salt water relationships exists. (6) In particular, contributions by Dutch engi- 
“neers have been most significant toward understanding these complex ground 
water hydrology problems. (4) Lindenbergh(5) has described the occurrence 
of ground water from that of the dunes ‘supplying ore 
) 


its history ry and problems. Detailed investigations 
Edelman (2) have dealt with evaluating the ground water distribution in the _ 

dunes in relation to planned ground d water ' withdrawals. A paper by Hui Huisman(3 

: ‘describes methods for determining permeability coefficients for the various — 

S dune strata. These coefficients have been adopted for the e analysis presented 
here and are included ir in the notation on the next page. 
a The following mathematical analysis enables the ground water movement 


i 7 and the boundaries of the fresh-water lens and pocket to be determined. “ies m 


dane and conditions too obtain solutions have been fit 
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reclaimed asa polder, the ground water level lies below mean sea level so 
that a landward flow of salt water exists. This introduces mathematical | 


_ complications which can only be solved with the aid of electronic computers. x i 


ot. Results presented below are based upon data for the area under study. vd 
_ Information on dimensions of the dunes, geologic strata, and permeabilities — 
was 3 available from exploratory measurements a and previous studies. . Although 
- the ‘quantitative results necessarily are not generally applicable, the method - 
of analysis is. Any coastal ground water situation involving a dynamic equili- 
= water over a stationary Salt water | body can be treated ina 2 
‘similar manner. In fact, with only r minor modifications for radial - flow, the 5 
procedure can be applied to analyze ground water — on an unconfined — 
The # authors have tried to emphasize the analytic ao involved by 
_ carrying a specific problem through to a quantitative end. It is hoped that the 


procedure will prove helpful to engineers who may encounter or relat-— 


ed -round water situations. _ 


OA cross- area of th the lens above 
Cy, Co, Cz, C5; = — constants 


2 average vertical thickness of the fresh-water name 


fresh-water excess of hg landward of the 


= depth of salt water above the 

= permeability of | the clay layer = 0.002 m/day 
= permeability of: the Holocene’ ‘sand formation, /m/day 


of th the ‘sand formation 

32 m/da 

: half-width of the fresh-water lens perpendicular to ‘= 


= half- ‘width of the fresh-water Pare including © 


N 


els = rate of effective precipitation, or 3/day 


D. 


fresh-water flow in the fresh-water lens, 


fresh-water flow in Pocket, m3/ 


= pumpage one side of the pe at x = 


| 
= g 
| 
= 

— 

— 
— 
a 

4 
water 
‘ound 
ment — 
d to fit ==. 


= horizontal coordinate measured from the 
the dunes perpendicular to the coast, m ed 

fresh- water ¢ depth in the pocket, 
= value of y y at the edge of the fresh- -water lens, m oe. 


= maximum value of ¥ at the center of the fresh-water 


er height in the from O O.D.. 


constant 


constant 
rat = constant =4 


=  piezometric level ( or of the fresh-wat -water 


= per of salt water = = 


average value of ¢, m_ >. 


‘pie. _ Although idealized, the sketch is a good approximation ee actu- 
al dunes cross-section before reclamation of the Haarlem Lake. The poeinons 
will be treated in two dimensions with flows referring to a unit thickness of _ 
one meter. The fresh-water lens lies above the clay layer, while the fresh- | 
water pocket with its een tongues lies below. All interfaces are boundaries 


.15 m O.D. (Amsterdam ordnance datum). "Salt water occupies ‘the 
above the clay layer on the seaward (left) side of the lens; however, on the = 
landward (right) side there is fresh water. Ghyben- ‘relationship 
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_ /.220-1.000 73) 0.2 26 
o 


“Therefore, a a constant V g gives a fresh- water elevation c of -0.15 + 0.2 baa 0. u 
_m O.D. on the landward side of the lens as shown in pe 


Analysis o of the Tongue 
“The tongues are those portions of the fresh-water pocket extending beyond a. 


on the landward and seaward side of the lens are symmetrical. Flow within a 

tongue consists of horizontal movement from the central portion of the pocket = 
an upward vertical movement through the clay layer. 
pu... the notation and Fig. 2, the tongue analysis can be manne on the ee: 
equation for flow within the tongue, it~ we, 


fresh-water lens. Because the salt water potential is constant, the — 


= 


POTENTIAL 


SKETCH OF ‘Dunes. cRoss- SECTION 


HY7 § HY7T 
the where pf and ps are the Iresn and sa ater densities respectivery. isthe q 
water depth above the cl 
than hg. Solving for Ah ang 
a 
fy 
and fee. 
at 
JQ __ Ke (F 0. 4) x 
At salt water interface the potentials in the two am 
But the salt water potenti a constant equal to -0.15 m, hence Eq. (5) 
— 
is of 720 = TO E = 


9 
Combining Eqs. (3), (4), and (8) yields the tongue flow equation _ 


+h "= 


wil 


Substituting 


With the permeabilities the of the can be ify, is 


known. This value will be obtained ‘subsequently in the analysis. pecuernt > a 


oY _ The form of the tongue interface can be ae directly from ‘Eq. (12) t by 


inserting values for Ko a and Kp; 


+63 = 


from ©, are also listed. 
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Aa. 


art 


— 


— 
(6) 
ae ‘ aa Cy is an integration constant. From the boundary condition at the end 
— 
— 
— 
oak 
— 
a »f the tongue, x = L and 
— 
— — 
g 
oJ — Ea. (12) defines the interfacial boundary of 
— = tongue L- { can be given by the boundary cond | i 
an 
«a es the coordinates of the interface listed 
x the obtained 1 
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4 Table 1 - Coordinstes o of the Tongue Interface 


= 
L-x 


37. 


101.0 


ae 


‘The fresh- ‘lens and are about the line of 
De - the dunes as shown in Fig. 2. _ Water is supplied to the lens by vertical perco- “! 
. lation from effective precipitation. , This equals gross precipitation minus eS 
J] cvapotranspiration and amounts to a rate N, assumed constant. Water leaves 
the lens either by outward horizontal flow at the edge of the lens or by verti- 
cal downward flow ‘through the clay ag 
In Fig. 3 is shown a sketch of half of the fresh-water lens cross- reection. = 
The downward flow the clay in the segment dx can by 


continuity equation for the same segment is 


iA 


Although D is a function of z, as a good it assumed 
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, then Eq. (22) 
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gs wiles Ce 
where C3 and C4 are integration constants. 
_ From the boundary condition that ees 


With this function of Zz, itis: not ; possible to : solve for y = f(x). 

necessary to find a good approximation of z = f(x) which al lead to a solution 

for y = f(x). Assume a relationship of the 


4 
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For’simplicity let A 
23) 
(23) 
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@ varies only a small 
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a or, upon integrating» 


center of the lens x= 


as 

or, , after integrating 


variable Q can be by y upon substituting in (3) ont ©), lead- 


2. o7 kp 


can be Four 
cases, different values of be considered and are listed listed in q 


6.3 


150 _ 


Substituting fc for from Eq. (6) Eq. (28). and inserting the second and 

(84) 


OM = 0.02 


HY 7 
ff 
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(30) 
| = 0, so that the integration constant a 
a 
- relat ater pocket inter- 
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values for x and from the third boundary condition (39), 


‘The ‘total from ‘the into the amounts to One half 


this quantity, or Qh, ‘flows through the clay layer between the center and the 7 
edge of the lens. Thus, from Eq. (31) with x = wi = = 3000 meters and _ 


0. 002 02 m/day 


7" 


where Q, is the flow through the section Yc into t the tongue. From 
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party Be 
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values of Q) Eqs. 46. (99) and (41), ‘respectively, 


- z= containing three unknowns. With these solved the functions y y= = f(x) 


By 

Eqs. (35), (37), and (42) contain the unknowns a, 7, Yc, and Cg. , The 
Cg can be eliminated by assuming a value for Ym> the four cases 
Table 2 will be considered. ‘This reduces the problem to the solution of three F 


and z = f(x) can be obtained directly from Eqs. (6), (28), and _. Results nei 
' the four cases are summarized in Table 3 and Figs. 4- ge i inte 


Table 3 - Values for “, 


0.092560 x 


|. 


With z= f(x) ‘now ¢ from the assumed Eq. it is to 


aie au compare this result with the more exact function according to Eq. (27). _ The F ma 


11 at x = 3000, so that from Eq. (27) 
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two unknowns C3 and N require two boundary conditions. The first of these dif 
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The second conaitio 


Eon. 27 
where A A isiaaiin the cross- -sectional area of the lens above O. .D. “Although the 7 
: integral from Eq. (28) cannot be obtained directly because ¢ = f(y) leads to an “ 
elliptic integral, it can be evaluated by planimetry. Thus, integrating ee re- 


| 


Eqs. (43) and (47) can can now be solved for the unknowns C3 and N. quantity 
oa appearing in both equations can be found from Eq. (45), as itsleft side 
equals No, (see Eq. (24)). The same is true for z, in Eq. (19) as a part of A. ae 
_ The resulting z-curves from Eqs. (27) and (28) are plotted in Fig. 8. A aia’ 
"maximum difference of 14 per cent occurs along the Case 1 curve; the — 
“difference for the four cases is 3 per cent. 
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Values of N for each of the — cases were obtained in the above analysis 


a a z. Each N indicates the average recharge rate from | precipitation minus» 
See necessary to maintain the lens and pages configurations. 


(Of the recharged a quantity Q percolates downward through the 
_ layer on each half of the lens, while the remainder N- -QQ flows outward at 
each edge of the lens above the clay layer. Values of Qh, found from Eq. (39), 
are listed in Table 4 together with N-Q) quantities. 
a = of N, Qy, anc and N- AN-Qg fe for the four cases are 2 also shown on Figs. 4-7. 
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_ Pumpage from the Fresh- Water Pocket 


q 
~ ‘If any water is extracted from the fresh-water pocket, the position ofits — 


3 4 interface will of course change. Under conditions in Netherlands water is “4 
si = at the edges of the lens (dunes) to minimize flow of water into the ee 


tongues where it is lost from a water supply standpoint. - The poepeen,: yn of th 
pocket interface can be determined for various pumpage rates atx=. a 


4 & _let AQ equal the amount of water pumped from one side of the ees at 


» Ye. any then be rewritten 
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| 
aa | 
_ ak. ‘2 
‘ 
1. (39), — 
3. 4-7, 
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80 that AQ m qs. (39) and (41) for a given value of y, Vo 
: 3 one of | the four previous cases. The value of Q) in Eq. (39) depends upon “ 
" and Y, which can be obtained by solving Eqs. (35) and (37) simultaneously. — 4 
Knowing | a, B, and Ce, the shape of the pocket interface can be calculated Bi 
from Eq. (33). Values of 4Q for selected values of E Yc} in the the four cases are 
and The length of the tongue with pumping is § given, as before, by knowing y,. Yo 
substituting in shape of the tongue interface is given by 


Positions of the pocket and tongue interfaces for values of Yo in Table 5 ? 
are plotted in Figs. 9-12. It can be observed that as the amount of ground ~ 
_ water AQ extracted from one side of the pocket is increased, the tongue depth 
a 4% _ ‘Ye decreases. The limiting situation of Yc = 0 in each case produces a maxi- 
; Bea mum water yield, which comes entirely from the pocket as the tongue no long- 
er exists. For ye =( 0, 0, hence AQ. By comparing values of 
4 naa in Table 4 with corresponding values of AQ at Yo = - 0 in Table 5, a substantial 


ow increase in flow through the pocket can be noted with the extraction. Quali- 


The at above analysis enabled the fresh-water flow in the ae and ein, as) 
3 well as the interfaces to be determined for specified boundary conditions. _ 
ae” ‘Today, drainage ditches in the dunes create a much more complicated ground 
, a water table than described above with z = f(x). , Further calculations are} 
ee. necessary to take into account the discrepancies. ‘Also the quantity N in actu- 
; od ality is more involved than described heretofore. It represents the difference 
5 
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causing an increased downward gradient between the lens and the pocket. 
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FIG. 10 “FRESH- WATER PUMPAGE FROM THE 
Rhine River water and water discharged from the dunes by evapotranspiration | 
and drainage ditches. There is no surface runoff fromthe dunes. 
_ Under peesent- day conditions the inland water table lies 5 meters below A 
_ mean sea level. _ This induces a landward flow of salt water under the dunes - 
< which shifts the interface to destroy the symmetry of the pocket and tongues 
F shown in Fig. 2. This factor increases the complexity of the problem to eof 
_ an extent that no direct mathematical solutions can be obtained. The matter = 


is being studied by engineers of the Municipal Water Works of —. 
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FIG. FRESH-WATER DISTRIBUTION WITH PUMPAGE FROM THE POCKET, 


Ir. C. Biemond, Director of the Municipal Water Works of Amsterdam, | 


q “kindly assisted in the planning : and general conduct ¢ of this s study. D. 


; the calculations required for solving the various simultaneous equations. 
Participation i in this study by the senior author was made | by a 
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Toronto, Assoc. Intl. d’Hydrologie Scientifi ue, vol. 2, pp. 168-182, 1957. 
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‘Lin | water from ad dune a Jour. Amer. Water 


S 


‘Todd, D. K. ., An of. to sea 
its control, Tech. Bull. 10, ‘San. Engrng. Res. ‘Proj., 
California, Berkeley, 72 pp., 1953. 


~~~. 


of 


| 


— 
— 
me a 
8 


ete 


> 


@ 


of the American Society y of Civil 


FOR EXCESSIVE RAINFALL 


_ thus conceal much design information. By relating probability to recurrence 
_ interval instead of rainfall depth the distribution of recurrence interval ae 
rainfall amount is obtained. This 


ig is 


a a function of depth. (1) The mean recurrence interval in years for a — . 

_y is defined by the equation R(y) = n/m, where n is the length of record oe: 

3 "years and m is the number of occurrences in the record equalling or exceed- — &§ 

ing the depth An occurrence equalling or exceeding y is called a y~event. 

a ‘ E.g., in a 20 year record the second highest value has a mean recurrence — a 
interval of 20/2 or 10 years. . The mean recurrence interval is also some- 

- times called the return period although the latter has usually been applied ae 

_ where only the extreme occurrence from each year is employed. Both ~~ 

_ frequency (or mean recurrence interval) and return period have the eumnd a 

_ that they give an average time interval of recurrence. SO oe —— 

The averaging of the time interval between y-events has always ‘cians aa) 
rather artificial, for the time interval appears to be the significant variable bi 

i and averaging obscures much of the information it might | provide. ty The effect ey 

7 of averaging has often been recognized by stating that the mean recurrence al 


e interval gives only the average recurrence period and coma — — will 


Note: Discussion 0; open until December 1, 1959. To extend one 
ae written request must be filed with the Executive Secretary, ASCE, Paper 2083 is ; 

a of the copyrighted Journal of the Hydraulics Division, Proceedings of the - 
_ American Society of Civil Engineers, Vol. 85, No. HY 7, July, 1959. - 
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fie 
Although recurrence intervals are of primary interest in Irequency analy- 
| 
ay 
The usu practice in ODtaining design data ior high intensity rainfall over 
\ 
— 


‘mean of time interval distribution is 


FE 

often occur less than one mean recurrence interval _— in time. ie For ex- 
Pot a mean recurrence interval of 5 years for a y-event gives no measure - o 
of the risk that y-events will occur in two successive years. This seems to 4 
= a rather troublesome deficiency, for the purchaser of a drainage system 
would be little impressed by i a mean recurrence e interval if the Sow of 


2 


curred in the first two years of operation. This deficiency will be overcome 
& defining a a probability distribution on the time interval between y-events _ 


Poisson Ditton 


‘ 


not 


e, 18 


probability P, ‘anda “rainfall less than y y with probability 1-p. “The trials are 

_ assumed to be stationary i.e. p is constant for all trials. The probability of 
: x y-events may be expressed by a Bernouilli or binomial distribution B(x;n ») it 
the usual x is the number of n is the total number 


each. If nis large, either an interval contains a Ly- -event or it does D 
not; hence if successive unit intervals are non-overlapping and independent, § 
the number of y-events x found in the unit intervals will be distributed in a ca g ft 
Bernouilli distribution B(x; sn,p). - Replacing the unit time interval by an inter- . 


fi a Bernouilli distribution with the number of y-events and non-y-events in each” 
t equal to nt to the nearest integer. The distribution then becomes B(x;nt,p). _ 
as _ It is well known that the mean of the Binomial vara teatnaiies is a= Leeds so the 


me small p, ate interest here, the Bernouilli distribution approaches the — 
Poisson distribution. In the case of the time distribution ¢ this may be aed as 
follows: Since (at)/nt, for small p, 1 - 


4 


For large = n the series approaches -at, 


Ry. 

— 

— 

Although most are familiar with the ordinary Poisson distribution, 
fa a Poisson distribution with a time variable, which will be employed herg 

as well known and requires some further explanation. (2) 
—— cae Consider a series of Bernouilli trials, viz. a set of independent trials with ee 

a 

4 1% 
4 
4 
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wie 


nm This is ‘the Po Poisson probability nai finding x y-events in a time interval t. Ese) 
‘The Poisson distribution of Eq. (3) is best fitted to observational data by 
~ estimating | a for a unit time interval. The nature of excessive rainfall and = 
its use in design data make it appropriate to choose t equal to one year. With 
» this assumption the parameter a has an optimal estimate a which is the ee 
arithmetic average of the number of y-events per year. 
The a~values for one hour duration and various y-events for Omaha, 
st a Davenport, and Minneapolis are shown plotted in the figure and tabulated in 
Table 4. These a- were in fitting the Poisson distributions 


| 
nter-— 
gives 


5 and one for t in Eq. (3). In ra tables x is the number of occurrences ofa 
tiers in the year, g, is the observed frequency, and P, is the observed 
probability. P, is the estimated from Eq. (3), &c is the 


“function. To examine the goodness of fit of the Poisson. distributions, Eo is 

compared with Simple inspection shows that there is reasonable a 

_ agreement, and this is verified by x2 values shown at the bottom of each © 
table. , All of these are exceeded in random sampling with a high probability “e 
indicating non-significant (N.S.) departure of the actual from theoretical 
frequencies and therefore a good fit by P,. Hence the Poisson distribution Rs 
provides a good fit to y- y-events of excessive rainfall. _ This is not surprising, 

_ for the excessive rainfall series quite naturally meet the basic assumptions _ 
for the Poisson distribution. For direct use P, provides estimates of the es = 
probabilities. for 1 various x-values and F. provides estimates of the Probability 
that x is equal to or less than a given value. These results, however, are 5 
only supplementary to the main objective— —the time interval distribution. pS 


} The Time Interval D Distribution 
ES Sou time interval t between : y-events according to previous discussion i 
= the random variable and has a distribution. This distribution can be foun 
) directly, but it suits | the ¢ mbectins here best to > develop it from the Poisson 
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abt. a From Eq. “) for the Poisson distribution 1 it is seen that the probability. 

For the distribution of t the probability is required that there is no y-event in 
the n intervals of length 5t comprising t and a y-event in the (n+1)th interval. 


This is is a joint probability given by the product « of probabilities of these events” 7 


which is the Poisson time interval distribution. The iets of this distri- _ 
- bution over ranges of t gives probabilities of time intervals between y- -events. 
Thus, the probability that the time interval shall be one year or less between 


ye ~events is is by the integral of Eq. a 


FQ) = | 
and in general time interval Tb by 


. ‘it may may be noted that only | the integral values of T equal to or or greater than one | 
§ need to be considered in the excessive rainfall problem. ‘Performing a 
By using a as the estimate for a, Eq. (7) provides all eeeianens coniee for 
recurrence intervals of length T. It should be noted that since the excessive — ‘zs 
rainfall series is independent and stationary, the good fit of the Poisson bests 3 
distribution to the data of Tables ‘ 2, and 3 also insures a a good fit to the time - 
interval distribution and no direct test of fit of the latter is necessary. | 
_ It is of interest to see how the time interval distribution relates to the “= Bi ie 
mean recurrence interval. Since Eq. (4) : is the distribution of the recurrence aa 
interval, its mean value should be the mean recurrence interval - Writing the -& 


usual equation fo for mean or r expected value value gives ia 


= 


discussion it was shown that a y-events occur on the aver-— 
«age per unit of time. If the time t is divided into a large number n of small _ _ 
— 
<= 
4 
a 
5 
= 
a 


= 


_the number of years divided by the number of y-events or ais whereas a be 


the number of y- y-events divided by the number of years of record or m/n. a 


ad 


Hence, a and the mean recurrence interval are reciprocals of each other as 

verified by Eq. (8). From this equation it is < also clear that much information 
has been ignored by only ‘computing the mean recurrence interval 1/a. _ The . 
a time interval distribution may thus be viewed as a tool for analyzing excessive 
rainfall data sufficient in itself, or as a means of extending an analysis which 

_ gives only the “mean recurrence interval. For example, Yarnell (1) gives mean 

- recurrence intervals called frequencies. They are estimates of 1/a inthe _ 
time interval distribution. Hence, Yarnell’s results may be extended to re- 


currence interval probabilities by ‘substituting the reciprocals of his frequen- 
3 Another aspect | of the relationship between | the mean recurrence interval 
and the time interval distribution is apparent in Eq. (7). Rearranging (7), - 
ing logarithms, and solving fora gives 


a i This setetes the reciprocal of the mean recurrence interval to KT) the proba- 


7. - bility of a recurrence interval of T or less and furnishes a means of — 


~ 


ing the mean recurrence interval in relation toa given F(T). = 


- a simi its properties it appears that the mean recurrence interval i 

a fairly arbitrary choice of statistic for design purposes. An event whichis § 
expected to occur with an average time interval of a certain length 
little light on a problem in which the main concern is for the individual time § 
= and particularly for the shorter intervals. . AS a consequence it ap- § 
pears that recurrence intervals other than the mean might be much more im- 5 
portant | than ‘was apparent in the past. Perhaps much of the popularity of the 2 
mean recurrence interval is due to its long usage rather than its real merit. J 

_ From a practical design standpoint and also from the view of the purchaser) 
x certain tain types of designed systems, the ‘shorter recurrence intervals might 


be the most important. “Where risk of annual recurrence is particularly im- 
portant the customer will be most anxious to be informed about this risk. In | 
_ view of this it will be of interest to give certain statistics for the one and two 

_ For this purpose the a’s of Table 4 were plotted in the figure and empirical 

ey. straight lines drawn to the | points on the assumption that the mean frequency — ¥ 
these ne exponentially. The smoothed a-values of Table 5 were read from | 
_ these lines for the various durations. nS The mean recurrence intervals which - 

= the reciprocals of the a’s are also shown there in parentheses. These — 

- agree well with Yarnell’s results(1) estimated from quite a different period of 


‘B) Using Eq. (7) with T= = 1 Table 6 was computed giving the probabilities of 
the recurrence interval being one year or less. For the lower rainfall depths 
oi it is seen that the chances of annual recurrence is very high, and itis not 


; until depths | of 2.5 inches are ventnt that the probability of recurrence is re 
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+The foremost objective of this analysis, however, was not Pence to refine 
“the mean recurrence interval or frequency analysis, but to present a method 
7 of analysis which might deserve consideration on its own merits. Forthis | 
purpose the approach is the natural one of all design problems, i. e.. , to set the 
_ risk or probability in advance and then determine the numerical value of the 
design variable. This is done easily by means of Eq. (9) and the figure. ~— 
Much more practical experience with the application of the Poisson time 
- interval distribution-must become available before its usefulness can be ap- 
_ praised. However, as illustrations of possible application two examples for 
the lowest two recurrence intervals, one and two years, ' which would appear — 
_ to be the most critical in design problems were chosen. In both examples 
_ 0.05 probability or risk of occurrence was assumed. Substituting 0.05 for 
F(T) and 1 and 2 respectively for T, a’s of 0.0513 and 0.0256 were obtained 
- from Eq. (9). These correspond to recurrence intervals of 19.5 and 39. 07 7 
"years. Entering the lines of the figure with these a values the estimated © 
| o ee _ depths of Table 7 were obtained. These are depths such that their one and 
re 3 two year recurrence intervals will only occur with probability 0.05. It is a? 
clear that the depths for two years must be less than those for one year since 
the former is a more stringent design condition. Pee ee 
oli From Table 6 it is clear that to insure with reasonable probability that the 
‘Spacing of values exceeding design (y-events) shall not be less than one years 


S 4 If the spacings between y- y~events are | to be less than one and two years a 
_ spectively with a probability of 0.05 or chance of one in 20, Table 7 gives the 
design depths which must be adopted. 5 It would appear from these examples” 
that in certain design problems \ where the most serious criticism arises from 
_ the too close spacing of failures, such as exceeding the capacity of storm _ 

a sewers, the only concern n might best k be for the probability of a cane re- 


y of the Des Moines Weather Bureau Office, for - compiling the excessive pre- § 

cipitation data, to Mr. Maurice Kasinoff for expert assistance in making all “a 
__, computations, and to Mr. Warren W. Buck, Jr., for preparing the Ts 
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normal distribution of is described. It is found that 
graphical procedure is much simpler and it produces results with less effort a 
than the analytical method. The method sm a herewith is illustrated ited with — 


—" the design of a hydrologic engineering project it is almost always found — =. 
necessary to know something of the probability of occurrence of hydrologic 
» events of various magnitudes for which the structures of the project a : 
be designed. From the viewpoint of safety, a structure in the project sh should By, 
be able to withstand the worst possible hydrologic event which might occur. = 
From the point of economy, on the other hand, the cost of the structure must | of 
be justified by the returns to be expected from its construction. . Consequently, - 
the selection of a design hydrologic event involves a calculated risk (or proba- 
bility) that a more severe magnitude of the event will not occur during the 
life (or design period) of the structure. In hydrologic design studies the __ ‘ted 
knowledge of the probability of occurrence of hydrologic events is investigated _ 
by a frequency analysis through the use of statistical methods. eave ee 
_ There are many methods of hydrologic frequency analysis in which the re- | 
corded magnitudes of a hydrologic event are taken as statistical variables and 
a statistical distribution of the variables is made to determine the “a 
characteristics, including of a given magnitude of the hydrologic 
event. According to Chow,(1,2,3,4) most methods of hydrologic frequency __ 
Note: Discussion open until December 1, 1959. To extend the closing date one month, ae 
se a written request must be filed with the Executive Secretary, ASCE. Paper 2084 is 
vr part of the copyrighted Journal of the Hydraulics Division, Proceedings of the me ae 
at American Society of ( of Civil Engineers, Vol. a? No. HY 7, ‘July, 1959. ; 
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and the frequency factor ofx for a given on probability of occurrence. 


a tie factor depends upon the ton of the statistical distribution of x. a i 
_ By elementary statistical procedures, the values of the arithmetic mean 
and the coefficient of variation can be « computed from the historical data of the 
7 hydrologic event. The determination of the frequency factor, however, will 
_ depend further on the law of distribution of the hydrologic event and the as- 
‘Signed frequency (or probability) for the design. With all these values de- 

termined, the magnitude of the hydrologic « event can be easily computed by (1), 

‘Conversely, if the magnitude of the hydrologic event is known, the frequency 

of the event can also be determined by (1). In this paper, the law of log- 

normal distribution is considered and the objective is to compute the corre- 

bution of the magnitude x is said to be log- -normal. The mathematical law 

which governs the log-normal distribution is called the log-probability law. 

_ It is generally found that hydrologic data of many kinds are log-normally dis- 

tributed. _A comprehensive discussion of the log-probability law and its 
_ heering applications has been given by Chow. (5) The determination of the 

eee frequency factor when a log-normal distribution is assumed has 


| 


oe 


_ been presented, (6) in which the method employed is purely analytical. | From ; 


these studies(5, 6) a number of { mathematical formulas have | been developed, — 
_ showing the relations between various s characteristic values of the variate x 


‘that is log-normally distributed. The formulas are found useful in 


bot 


in which x = the magnitude of a hydrologic 
— 
a 
| 
iz 
a 
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a 


FREQUENCY FACTOR 
s the additional notation is as follows: | 
the coefficient of skew of the variate x 


= logex, | the natural of x 


= the x whi which oc occurs -of - time, 


‘The pr problem is to determine the frequency factors of 


the log-normally distributed variate x at various assigned probabilities when | 


the value of either C, or Cy is given. The graphical procedure 2 this de- 4 7 

5: ® When C. is given, compute Cy by trial and error, using (2). 4 It can el my 


; (2), when Cy is given. 


Compute M by (7) x = 


(4) Compute xt by (9), 
(5) It has been shown(5) that the probability at mean is that which cccurs q 
when the frequency of x is equal te to and this is 
Jr) 


y=9 y/2, the of the probability Find the > 
value of Py, which may be obtained fi from a normal function 
table for Ky = oy/2. 
Ona probability paper such as shown in Fig. ‘represent 
_ the ordinate by x/x, in which x is taken as unity, and plot any two of the 
_ following three points: (Ordinate = M, abscissa = wo (1, Py), and cl 
rs (7) Draw a straight line passing through the plotted | points. = Hip Shi 
On the plotted straight line read the value of x/x corresponding to 
a) For known values of Cy and x/x, compute the corresponding frequency — gq 


illustrate the above procedure by a numerical example, assume 
Cy = 1.139. Then, by (2), Cy = 0.364. By (4), cy = 0.3527. By (7), M = 0.940. 
By (9), x* = 1.338. Since Ky = = 0.3527/2 = 0. 1764, the normal probability as ie 


,  Drawa straight | line on the log-normal probability paper in Fig. 1, passing oe 
; through the points: (x/x = 0.940, P = 50%), (1.000, 43.0%), and (1. a 


» 15.87%). Actually, any two of the three points are sufficient for the con- 
struction n of the st straight line. line. | 
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“for equal to or greater than the value of 
‘Using with Cy = 0.364, the factors are computed a listed in i 


Table Ps of of frequency factors” for 
= 1.139 and C_ = 0.364. 


2 y im percentage oo 
to or greater than the given 
x) 0 00414 0. 0.700 0.3 ‘sank 620 2 1135 209 
x) -0.684 -0.586 -0. 0.472 -0.300 0-265 0.68 0.680 1. 01.135, 1.797 2 

= -1.879 #1.297 -0.824 -0.165 0.728 1.869 3.119 4.997 6.624 

2 A table of approximate -normal frequency factors was first prepared 
Y : by Hazen(7) and then revised. (8) Until recently, a corresponding table of ae 
exact log-normal frequency factors was computed by an analyti- 
cal method. 5) It is believed that a great amount of time and effort can be 
saved for either preparing or expanding the table by using the 
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L WAVES AND SLUG | ‘FLOWS! IN INCLINED OPEN CHANNELS# 


‘Paul G. Mayer, M. ASCE 


1 


pie and are studied phenomenologically in an inclined open channel. The 
basic characteristics of flow are expressed in terms of pertinent physical es! 
properties. A theory is presented and reference is made to similar phe- _ 
nomena in allied fields in which mathematical analyses exist. It is hoped that 
a study will — the present knowledge of unsteady a in open , 


ie in open n channels is characterized by a free ial with constant 
"pressure conditions. Possible steady state regimes have been designated 
conveniently as laminar- ‘subcritical, turbulent- 
subcritical and turbulent- -supercritical. However, unsteady flows exist 
are not — to analysis by methods applicable to any of the — re- 


One such regime occurs in 1 shallow flows in inclined 
channels. It is characterized by intermittent surges and wavy patterns 
has been variously called roll waves, rain waves, slug flows, and other names. of 
Hydraulic engineers encounter these unsteady flows in inclined open channels © a 
and spillways. The increased height of the waves requires additional free- 4 
board to prevent spillage. The concentrated mass of these surges calls for vie 
added structural safety factors. against transient pressures and stresses. In ey 
Note: : Discussion open until December 1 1959. To extend the closing date one month, | 


a written request must be filed with the Executive Secretary, ASCE. Paper 2085 ~ = 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of 


the 
American Society of Civil Engineers, Vol. 85, _ HY 7, July, Sonar 


a. Presented at the August 1958 ASCE Hydraulics Division Conference in 
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— chemical engineers, the phenomena are important in liquid- gas reaction 
processes. Liquid mass transfer rate is basic to diffusion reactions and the — 
rate is greatly increased in unsteady 
As a result of this investigation, two distinctly different wave ‘train phe- a 
nomena can be delineated. They are designated at a as roll waves” 


forces. at 2 a slightly disturbed surface. ‘They a are . characterized by transverse 
ridges of high vorticity and intermittent quiescent zones (Plate I). 
- Slug flows result from instabilities which cause the transition from super- 
critical laminar to turbulent flow. Locally « disturbed regions spread trans-— 
versely and contaminated adjacent zones similar to turbulent spots in wind a 
_ tunnel investigations. | Slug flows are characterized by a succession of highly 
‘The Theory | of Roll Waves | and SI 


to both surface and gravity forces. wong 


own When in capillary and gravity waves result. The general 


Sy equation for the celerity of surface waves contains expressions for both 
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- the wave speed relative to the medium. A minimum celerity exists for surface 
waves when both influences are equal. Cmin is obtained by differentiation of are 
the General wave in respect to A. Then, ‘ 


_ When small disturbances are impressed upon a surface, the celerity of each 
"wave is is determined uniquely : for every wave length. ‘However, wave groups a 
_ are known to travel at speeds different from the celerity of solitary waves. 

group | velocity of capillary waves is 


and for waves 


and fous are and do not conform strictly 
ee. In open channel flow with a disturbed surface, both capillary waves and - 
gravity waves exist. Their behavior and the competition between them are re- ar 
; sponsible for the formation of roll waves and slug flows. _ Most solutions of ce 
instability problems and the analysis of progressive waves are based upon the 


assumption of parallel flow. 3 In this study, the proximity of the boundary is be fs i 


important and the initially | smooth flow suggests viscous flow. For a case of 


one-dimensional viscous motion, the Stokes reduces to ty 


ay 


2 


#) 


a ves | Ss 
and in shallow flows, ‘their approaches y gD. Celerity is defined 


iii 
_ Upon integration, the velocity becomes 


slop opes 


For progressive waves to form, the following 


‘For con ts wh equalit 


— 
— 
atisied 
viding by y/gD, whe 
4 Interestingly, this critical value number was also established by 
 M.T. Lighthill and G. B. Whitham,‘3) who used the Chezy equation intheir 
i r _ Other significant parameters are the slope and Reynolds number. Using f— 
velocity and d s numbers become J 
ies = critical Froude number can now be interpreted in terms of a Reynolds — ' ai 
number and a slope. From expressions (12a), (14), and (15) 


Squaring both sides and simplifying, ‘the relationship becomes 


and in in terms of Reynolds number and slope 


VR, 
which leads at NF = 2 


g ‘ 
As a interest, the formula by Dre 


s 
where C, was the of Chezy’s Yo, ‘the 
be expressed for flow as 
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ia criterion was established by assuming the surface fluid 
ss velocity equal to the velocity of a small \ wave so that the profile remains un- 
_ If the Froude number is less than ate 2, the wave velocity exceeds the 


surface velocity of the stream and 


“creases in height. “Also, the higher seen of a wave profile move with neil 


«Fig. 3 


Rig. 4 shows how long waves decelerate as forces 


As the wave length decreases the wave ‘height increases and neeaking of “= = 
: progressive waves in shallow water always takes eae LL eg 
; 7 Perched _ The frontward steepening of gravity waves and the subsequent increase of 
capillary effects are primarily responsible for roll wave ro Two 


Long waves steepen, amplify and break. If the waves 
are influenced equally strong by capillary and gravity forces, a wave train " 


ing of the waves s takes p ang a spectrum of wavelets re- | 


. Each in turn is to surface and forces. ‘The | are 


phen 
— inusoidal pulses : 
turbances may be sinusoidal p 
— stur 
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assumed to be of the same height but not necessarily of the | same heenanen. = 
‘ae the wavelets may have celerities different from the initial wave, some 
are assumed to interfere with each other. This interference is similar to tl the 
a "phenomenon of sound waves which have different frequency and phase re- -— 
lations. The resultant of two waves of nearly the same frequency varies in 
amplitude from zero to twice the amplitude of either component. In the —_ 
"channel, the process of interference may repeat itself until waves oflarger _ 
fe amplitude and selected frequency predominate. Finally, the forces of — 
: exceed the surface tension effects. The larger waves overtake smaller waves 
= coalesce with them and grow even geen ‘After ; a formative period, roll 


Wheat the Froude number is larger than two (NF > 2), the velocity ofa 


small surface wave is less than the surface velocity of the fluid. The a. 
ing of the wave occurs at the upstream end. The ~egetrauaned in curvature causes oes 
a decrease of wave length and a deceleration of aun wave, and consequently, Fos 
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the wave velocity is less than velocity, the higher velocity 
fluid enters a region decelerated by the breaking of the wave. Consequently, 

the main stream is also decelerated. Some of the kinetic energy is converted 
3 into potential energy and the depth of flow inc reases. Some of the kinetic | -F 
_ is expended in . agitation and | mixing. _ The combined effects are a local 
_ hydraulic jump which marks the transition from sua reir flow to 


Boas a The breaking of the wave unbalances the dynamic equilibrium of the super- 


> 


critical stream which is potentially ready to become turbulent. In laminar 
supercritical flow on a smooth channel, the breaking of waves triggers the | 
a transition to turbulent flow. The change in flow conditions exhibits a sudden 
4 increase in resistance to flow. Although the channel slope and the boundary 
et are unaltered, the high velocity flow enters upon a friction slope 
not sufficient to maintain that velocity and a hydraulic jump occurs. 
Subsequently, the agitated region 1 grows into adjacent undisturbed flow. 
_ The velocity distribution in each regime is different and mixing takes place. 
Vorticity transport and turbulent agitation cause lateral contamination. tf ig 


4 
‘ 


stream (Fig. 7). When an earlier jump intercepts this high velocity — 
: _ Stream, the jump front remains | a negative bore until all available energy is 
f 2 fed into it. Then the jump decays an and forms an expansion on wave and the down- 
stream face of the disturbed region becomes a positive bore. Subsequent _ 
_ bores telescope and form the ‘Chatacteristic saw-tooth cou of - flows 


(Figs. 8 and 9). 
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counters great difficulties. Convenient approximations often lead to 
of questionable veracity. ‘The background of experimental reality provides — = 
important insight into physical phenomena. Interpretive physical reasoning 
and dimensional analysis are other tools important to aninvestigator. = 
The variables of a phenomenon may b be into 


a) Boundary conditions, 


initiation and subsequent development of roll waves es and san flows were: 
. depth at point of initiation of disturbance 
Dy depth of average wave crest 
4 channel station at point of ein 
slope of channel 
volume rate of flow per unit widt 
‘surface velocity of water — 
‘velocity of roll wave 


Red 
: Vel: » velocity of slug flow 


urface tension ane 


viscosity 


flow characteristics and boundary conditions were in the 


‘Fluid properties were taken fro from tables. 
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in n this. xt included the test flume; equipment to 
measure pion yea and velocity; a disturbance generator; and — 


The experimental work was carried out with a tiltable flume (Fig. 10 

Plate II). The channel consisted of two sections. Each section was 12 feet — 
E long, 18 inches wide, and 12 inches deep. The channel floor was covered yO 
P. two 1 12 foot long slabs of 1/4 inch thick plate glass. ‘The s side walls, also of : 
plate glass, were installed in 4 foot sections. The underside of the floor slabs — 4 
were painted black. The side walls were left clear. All longitudinal and ~ 5 ahh 
joints: were cast with paraffin and planed down to con- 


ditions. 


a 


Water entered the test flume t Miiensih a 3-1/2 by 2 by 1-1/2 foot head tank. 
Wire screens were placed in the tank which supported a graded gravel and  _— 
sand layer four inches thick. This filter course was intended to reduce the 
_ initial turbulence in the flow. A curved steel section extended into the stilling _ 
was aligned with the glass hottom in order to obtain smooth entrance 
"conditions. top members of the channel | frame served as supports 
 — brass rails of 3/4 inch diameter. This track facilitated the —_ 


one foot stations, at the crest. 
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The water supply was mensared by a ser 
area flow meters. Their "respective effective were: 
From 0.001 cfs/ft. to 0. 006 cfs /ft. 
x ‘b) From 0.005 cfs/ft. to 0.030 cfs/ft. 


¢) From 0.010 cfs/ft. to 0.060 
- _ These | meters had been calibrated and were checked intermittently by collect- 
ing the flow for | gravimetric ‘measurements. A Toledo | dial type 
was used. All water was subsequently wasted. 
Depth measurements under steady flow conditions were made with a point 


«wage. ‘This, gage ‘was mounted on the instrument carriage. For greater 


a fay, an attached Federal dial indicator allowed readings to 0.001 inches. 


For depth measurements in wavy and unsteady flows, the variations in 
hei height were recorded by the electronic wave recorder. The primary com- ~ 
a _ ponent of this consisted of two 0.005 inch thick platinum probes separated a: 


a distance of about sixteen millimeters (5/8 inch). The passing water ented 


ie > ~ the electrolyte. The signal was amplified in the wave height measuring s 

: device. This amplifier had been built by the Instrumentation Branch of sal 
U.S . Army Corps of Waterways Experimentation Station, 


q variations in depth (Fig. 1). 


5 “The Brush oscillograph could be operated at three different ot enti 5 mm. “et 
sec., 25 mm./sec., and 12.5 cm./sec. The slowest speed was generally used 
to measure wave heights and frequencies. The faster ‘Speeds were more able 
me For the determination of surface velocities and wave velocities, the time yl 
_ measurements were accomplished with an electric clock capable of recording 
0.01 sec. The clock was activated manually by n means of a microswitch. _ The 
_ switch was attached to a long extension cord and thus permitted utmost mo- oa 
bility. The channel stationing was used for length measurements. cheery 
ae A small bottle of compressed gas was attached to the head wall of the still-_ 
ing basin. A flexible hose connected the bottle to a control valve located = ae 


could also be regulated by the valve attached to the gas container. 
_ All pictures were taken by the author with a Kodak Graphic View camera. 
Royal Pan 4 by inch sheet film used exclusively, 
_ The laboratory procedures evolved during a period of preliminary work. Te: 
precautions were detected early and followed in subsequent 


ments. Precursory investigations delineated the ranges of slope anddis- 


and flow characteristics \ were observed and recorded repeatedly. 


a Five different channel slopes were used. Their respective tangents 
— 0.0175, 0. 0349, 0. 0524, 0. 0699, and 0. 0875. On each slope, a series of nine 
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The oscillograph records were readily convert- 
ith the ain of coalihration curve The accuracy _ 
a 

3 
al 
| 
a ; _ After the channel slope was set and a given discharge obtained, the instrument E ae 
Pes a | carriage was moved from station to station. Considerable time was allowed _ a 
ne | to elapse after each movement of the carriage in order to eliminate external [i 
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as as possible. Depth measurements and wave 


_ were then made. Test results obtained under these conditions were r recorded . 
— mean velocity of the ‘Suaten eh water | was calculated from the known con- 
ditions of discharge and depth. . Surface velocities were measured by aos % 
im 
_ specks of talcum powder to travel between given stations and by measuring Pas 
elapsed time. This was repeated four or more times for each 


Wave velocities were established by 1 measuring the time interval in which 

a wave crest moved between predetermined stations. 

a ‘The entire sequence of measurements was repeated with the disturbance 

generator in Test results under these conditions were reported 
4 _ ‘The initial location, type, and subsequent development | of the surface dis- 
turbances were observed and recorded for every ‘setting of slope and dis- q = 


_ charge as well as Regime *Q” and Regime *S”. Photographs of significant = 


events were waves apes photographic effects were obtained indirect 


_ complished by a subdivision of the material into a logical sequence. Nj First, he 
_ the flow conditions leading to instability are discussed. ‘Then, roll wave and — 
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‘Unatabte to ‘one Formation 


aT, 
in the smooth laminar | sheet flow degenerated 
Pos ‘toh a flow marked by transverse ridges. The instability occurred in a flow _ 
which was initially | undisturbed. o can be attributed toa natural mode of = 
os 
_ the superposition of finite disturbances by means of the disturbance gener-_ na oe. 
B ator. _ The distinctive difference between the operating conditions was that the 
Bee - breakdown of the flow took place further upstream when finite disturbances 
_ _ External disturbances of various amplitudes and frequencies were used. 
i Both factors seemed to bear upon the stability phenomenon. “However, — 
_ breakdown occurred sooner with larger amplitudes. No quantitative measure-— 3. 
in of the disturbance amplitudes or frequencies were made. Significantly, 
the instability of flow, and hence the formation of roll waves, was eliminated ae 
_ when a wetting agent was added to the water. A wetting agent may affect not hts 
only the surface tension the liquid | ‘but also cohesion and spreading. Since 
ecsuies that of cohesion, an inter-relationship exists between surface . 
- cohesion, and spreading. (9) 17 The exact action of a specific wetting agent upon 
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_ mation of roll waves. The fact that instability of flow and formation of roll __ 
waves were not observed when a wetting agent was added to the fluid is helpful ‘a ; 


to hydraulicians engaged in model work. 
4 After the breakdown of the laminar flow was precipitated, the development 


caused the e instability of flow. Hence, laminar flow in inclined open achannels = 
was unstable because of inherent disturbances as well as 
perturbations. Surface tension played an important role. 


“Unstable Flow Leading to Slug Flow Formation 


‘aa as the transition from laminar to turbulent flow. Instantaneous point ; 
disturbances occurred randomly in the smooth sheet. disturbed regions 
spread transversely and were swept downstream at the same time. The a 
, photograph on Plate VI shows this phenomenon. The highly agitated regions — . 
differentiated the turbulent flow from the surrounding stream. The awe of fe 
_ these agitated zones was similar to oblique traveling hydraulic jumps. In the | 
_ literature of aerodynamics, these transitions are known as turbulent spots.(10) _ 
_ Superimposed disturbances from the disturbance generator had no appreci- ie 
“able ¢ effects upon the stability of the flow. However, raindrops | onto the = 
channel and sand grains in the channel precipitated turbulent spots. The 


Analysis of Experimental Data Initiation Initiation Conditions 


‘toasbeaman boundary and flow characteristics of the undisturbed flow ine 


Length Reynolds | Number, I 
Number, wy wind 


_ 


the presentation of the various characteristic curves, specific symbols — 
were selected to represent certain flow conditions. A legend =r the — 


ox _‘Tipple flow, train of waves with short wave length, 


roll waves, 


waves 3 at of breakdown, usually a at Rp ~ 420. 
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with the slope as the third parameter. An analogous relationship was obtained | 4 
d = when the gravitational and the surface forces were compared (Fig. 13). The g 

i transitional character of these curves is obvious. The graphs on Figs. , 12 and 
= 13 show that the surface tension and viscous effects predominated in the phe- | 

j - nomenon of roll wave formation. At higher Reynolds numbers, corresponding — 


to the conditions which led to slug flows, viscous and surface tension effects _ 


4 decreased and, finally, became negligible, 


= _ From the theoretical considerations, it was seen that the mean velocity 
& _ and the depth of flow were a function of the slope. Thus, a simplified parame-- 
a 5 ter, the modified Froude number, NF/ ys, and the | Reynolds number could be 


sentation of Froude number, Weber number, and slope (Fig. 15). a 


: _ A combination of all parameters was undertaken so that an overall grasp 
of the instability phenomenon might be obtained. | This also has practical engi- 
importance. In practical one m may know the conditions of E 


temperature. dimensionless parameter r containing the rate of flow, a, 


Was named “critical flow mumber", For the one containing the depth of flow, 
was named “critical depth ni number” » Np. ‘The derivation of these 


_ The instability condition was a function | of Rp Nw, and 
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plot of Fer and Np is in Fig. 16. In the region instability 
which led to the formation of roll waves, the two parameters showed a ee. 


Pa 


2 ‘The “critical flow number” simply states that instability is imminent when 
ever the surface tension forces become comparable in magnitude to the mo f 
- mentum carried by the stream. This is equivalent to saying that the free 
surface energy becomes comparable to intrinsic energy of the stream. a 


| The “critical depth number” is essentially a Reynolds number for laminar 


flow. The expression was developed earlier and given as Eq. (28). A eS 


| dias expresses the ratio of kinetic | , energy in the stream to the energy dissi- 
= by viscous shear. Thus, the For ant Np» have perfectly 


of 


‘ Fy parameters when a parabolic velocity distribution was assumed. All experi- 1 
g mental data were plotted in Fig. 17 to show the general inter-relationships. ea Tt 


ae The instabilities which led to the formation of roll waves were in the region - 


nis of supercritical laminar flow as predicted. Some instabilities leading to = 


aa flows satisfied the predicted conditions, but most data indicated that slug 
flows were a transitional phenomenon and the actual velocity. distribution may — 


tak, ne (17a) and (19) predicted that: ad 


Sey oe Roll. waves were observed up to Reynolds numbers of 410, 413, 419, and _ 2 
é 418.0 Then, for a mean RDmax = 415, the minimum slope for the formation 
Toll wa es was, according to (1 


ri 


— * 
| 
Be 
: a Experimental Relationship Between Reynolds Number, Froude Number, and | x 
if 
4 ar 
gr 
lal 
e 
4 
‘which was verified by the experimental results. The corresponding Froude ~ 
‘number was experimentally equal to NF = 2.06. At the higher Reynolds 
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matiiia the flow often degenerated into turbulent spots below the overflow a 
| entrance to the flume, before uniform flow conditions were obtained. Instabili- _ 
‘ties which led to the formation of slug flows were usually observed at Reynolds 
numbers larger than 1250. Turbulent spots were produced below this value by . 
artificial rainfall. . The deviation from the theoretical curve > (Fig. 18) also va 


starts at a en number of 1250. With a theoretical minimum Froude 


patel: slug flows were re observed experimentally on ona lope of — 


ocation o ace Ins ity 
ts It was observed that the location of the initial disturbances senting’ to ase” 


wave formation moved downstream with increasing Reynolds 


ae 
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, Reynolds ‘mabe and Slope 
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“number” and the “depth Reynolds number” elucidates this phenomenon 
g _-—-—s- The plot does not reflect the locations of instability for artificially dis- 
_turbed flows. Finite disturbances from the generator caused the instabilities 
to occur further upstream in the case of roll waves. For slug flows, the ef- 
fect of the generator was less pronounced. | ‘Whenever the Froude number was 
7 excess of the theoretical minimum of Np= 2, turbulent spots could be :, el 
es ik _ brought about by raindrops onto the channel or by roughening the channel bed 


ss prediction of the location of the initial surface disturbances cannot be 


Wit = 


of Roll \ Waves at) of 


= short wave length. The speed of these waves was always larger than the © 


ay: 
where the minimum , experimental conditions was: 


over the entire length of the test flume. 
q = - Roll waves did not form on . slopes below 3%. The formation of roll waves 

_ is shown schematically on Fig. 20. Perfectly smooth laminar sheet flow first 

maffered very slight undulations. breakdown of the flow into ripple flow 
™ followed immediately. Then, : a confused zone of rearrangement preceded the 
| aa _ final emergence of roll waves. Plate I shows a photograph of this process. — ei 
— waves were characterized by transverse ridges of high vorticity. The >? 

zones between the wave crests remained Iv, Vv.) 


= 
Th f Sl 
Formation of Flows 


Slug: flows originated with local disturbances. | 
= ‘transversely and were swept downstream at the same time (Plate _ A = 
schematic representation of spot growth is shown in Fig. 21. A i 
wea made between the experimental results of wind tunnel investigations ‘aa 


.- development of the disturbed regions in the inclined open channel. The 


 quvelopes of spot growth behaved similarly for Schubauer and Klebanoff’s(11)_ 
a tunnel measurements and Amein’ g(12) result in the test flume (Fig. 22). 
_ The regions between the spots were quiescent. _ The upstream front of the | 


disturbed regions moved with a velocity approximately equal to wt of the 


. stream front resembled a traveling oblique hydraulic jump; the downstream 
_ end developed into a bore. Slug flows resulted when the quiescent pockets" 
were eliminated by successive bores which telescoped » the ot 
charact eristic saw-tooth profile (Plate VoL.) 
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Plate VIL. History of Slug Flows 


he 
spot incidence increased as the flow progressed downstream. An increase in - 
spot frequency was also observed with increasing Reynolds numbers. The — 
disturbed regions grew as previously described. A drawdown of the water ie 
surface corresponded to a decrease in flow. This observation corresponded : i. ; 
to the “slackening” of the current reported by Cornish(13) halfa comny ago. 
= Slug flows were characterized by highly turbulent intermittent surges. 
Transverse ridges separated the regions of agitated flows and numerous _ 
smaller ‘waves were seen to ride the larger | surges. During the formative 


custo wave fronts. ‘Eventually, the waves into a common trans- — 


Wave Characte ristic 


; ty were investigated. A comparison of the data of wave height and wave ve- 
locity showed that a terminal roll wave velocity existed. The relationship ._ 
| plotted in Fig. 23. The terminal velocity \ was then related to channel a and — 


plot, the r resulting equation was: oso ae 
Terminal Ig= 
vie 
On the graph, the surface velocities were It was seen 


The orderlines suggested that a dimensional above was di might be ciate. In _ 
} the following analysis, both roll waves and slug flows were treated alike. ‘a 


mensionless parameters: 
te 
Reynolds — 


Froude n number, 

_ Dimensionless roll wave 


| 
‘Depth 


Wave number, 


The physical of wave number pony be “kinetic 
energy per wave per unit discharge”. 


usual symbols were employed. 


¥ 24 pointed out that the roll wave 
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4 Reynolds number to the one-third power. A plot of dimensionless velocities 

4 ae j and Reynolds number showed the same trend for roll w: — 

ia | lows, the velocities were always less than the surface .: — 
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in dimensionless form, showed ia constant r ratio independent ynolds 
: When all the experimental data was , plotted in terms of the modified wave 
_ number and the ratio of wave height to wave length, a narrow region delineated 
the zone of wave formation. The data pertaining to the steeper slopes plotted 
_at the lower portion of the loop, so that the characteristics of roll waves and 
slug flows on still steeper slopes can be fairly well predicted (Fig. 26) rk 
_ The plot involving the wave number pointed to the fact that this parameter 
tends to zero with increasing discharge. Slug flows were characteristic at _ 
oe discharges and plotted below wave numbers pertaining to roll waves. ‘ 
At very large discharges, the wave number vanishes. ‘This reflects the fact 
that roll waves and slug flows are phenomena observed only at relatively low 
discharges. Although the continued growth of wave height and wave length 
made the length of the experimental channel seem insufficient, the pies of 7 | 
bs height to' wave e length approached a limiting value (Fig. 26). 
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nstability were either carried by the stream or were ‘externally 


a Roll waves formed whenever the surface velocity oft ‘the ewer stream 
Pe was less than the minimum velocity ‘of surface wav es, provided that the ag = 
ie channel slope was sufficiently steep. Experimentally and analytically, ei 
a a minimum slope sufficient for roll wave formation was about 3%. The maxi- 4 
‘mum Reynolds number of the undisturbed stream was 420. Roll waves  result- 
_ ed from the interaction of surface tension and gravity forces at a slightly dis- 
turbed surface. They were characterized by tr transverse ridges of high vor- 
tieity. regions between the crests were quiescent. The terminal veloci- 
ties of roll waves were in excess of the surface velocity of the undisturbed oF 
- stream. . The ratio of wave om to wave length approached the limiting value 
- Slug flows resulted from instabilities which caused the > transition from “ 


a hydraulic jumps and seemed — to the turbulent spots - 


a were not considered to be “shock” waves. "cue one may oT 
Le shallow flow in an inclined open channel a boundary layer phenomenon, the | 


Be ® Slug flows formed whenever the surface velocity of the undisturbed flow 
was in excess of the minimum velocity of surface waves. They were charac- q 


tative analogy between compressible and incompressible fluid HOW. ow 


Te - Reynolds numbers of 4000, the resulting flow was | thoroughly turbulent and no 
slug flows were observed. The minimum slope for slug flow formation was h 
- approximately 2%. The Froude numbers were always in excessof2. The | 
flow decreased noticeably between waves and the surges literally jumped o off | 
the experimental flume. ‘9 In prototype channels, slug flows are known to occur 7 


with the same ‘speed. The ‘@eviations from ‘the m mean an were a at times ; considera 
ble and the concept of waves “that are alike in size and shape and travel with 
a constant velocity” ° which was expressed by Thomas , Dressler, and other “7 
_ writers must be questioned. In this study, the analysis of the oscillographs as 
well as the measurements of wave velocity reflect average values. 
__ Roll waves and slug flows were - analyzed and plotted in terms of pertinent 
‘dimensionless parameters. Flow conditions susceptible to instability were 
_ delineated by a “critical depth number” and a “critical flow number”. Signifi- 
cant wave characteristics were expressed by a “wave number”. 
—aaee boundary conditions enhanced the formation of roll waves. 5. Rain 
waves, the often observed on street eurfaces, are 
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waves, and it seems they will form, of the 

boundary conditions. _ Experimentally, slug flows were also precipitated by 
_ boundary roughnesses. ‘However, the disturbance sources were isolated three-_ 

_ dimensional obstructions (sand | grains). A . sufficiently large number of dis- 
turbance sources, that is, a New rough channel, agitates the flow sufficiently 


of swift mountain streams. The snubieatie of the channel as a measure to . 
eliminate slug flows was pointed out by H. Rouse.(14) 
| he _ The theory of roll waves and slug flows presented in this study presupposed a 
_ that the velocity distribution was significant. A parabolic distribution was __ 
assumed and was found quite satisfactory. Jeffreys(15) reported earlier - } 
_ slug flows (alias roll waves) were the transition from “turbulent flow with a a : 
_ nearly plane surface to turbulent flow marked by transverse ridges”. _ ae 
_ The pertinent characteristics of roll waves and slug flows were ee ae - 
in terms of dimensionless parameters. In order to test their validity, a simi- 
lar study could be undertaken with a different fluid. A longer c channel capable — 
of steeper slopes would also extend the range of this study. —es ha 
: * ‘Disturbances leading to the formation of roll waves could be be superimposed 
a the laminar stream by controlled m means. In the | experiments, randomly — i 
_ generated disturbances led to the same periodic wave formation with the same 
selected frequency that evolved when an undisturbed stream degenerated into ‘ee 
roll waves. The frequencies and amplitudes preferentially leading to roll 
| wave formation could be studied by means of a vibrating ribbon or controlled a 
electric pulses. ‘Similar ‘methods have been 1 employed ‘successfully in wind © — “ 
Slug flows were reported to lead to air entrainment on dam spillways. 
Maistre and R. Oblensky(16) published photographs which Suggested 
air entrainment was caused by surface roughness. Thus, it seems incongruous — - 
that rough channels eliminate slug flows but generate air entrainment, if air end 
entrainment is preceded by the former. A further study with a rough channel 


is planned and could elucidate also the proper significance of the two phe- - 
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by Martin E. Nelson and Johnson 


MARTIN E. NELSON,! F. ASCE, and HARVEY J. JOHNSON, 2 M. _ASCE.— 
nt Webster’s discussion showing a good correlation between model and a) << 
——- hawser forces is a valuable addition to the field of expe rimental Ba * 


with those in the basic paper. In some instances the hawser force per ton of 


tow or ship displacement can be compared aut in this instance the lifts are too My 


location o of tow in the lock, rates of filling or emptying, depth of water cushion { 
unter the vessels, and the elasticity ofthe restraining hawsers. = 
The senior author has observed the operation of the prototype Eisenhower 
and Snell (formerly Grass River) locks. The water surfaces in the chamber © 
and lower approach were very smooth and the valves functioned quietly. In a 
_these respects there was good agreement between the model and prototype. te 
_ Although accurate tests have not been made in the prototype locks, it appears | . 
_ that their hydraulic systems are somewhat more efficient than their respective 


_ models indicated, the time required to fill or empty the locks being somewhat 


shorter. ‘This type of discrepancy between model and prototype performance 
which has been observed in other locks and in outlet structures, is probably 
“due to lack of exact similitude in conduit roughness. The vortical activity a 

the intake manifolds was more in the full | scale locks" than in the 


models. 3 


+ 1. ee. Engr., St. Paul Dist. Corps of Engrs., U. S. Dept. of the Army, Rei 

2. Hydr. Engr., U. S. Dept. of Agriculture, Agriculture Research Service, ; 

q a Soil and Water Conservation Research Div., , St. Anthony Falls Hydr. Lab., B 

_ Minneapolis, Minn., the Corps of Engrs., Ss. of the 
Army, St. Paul, Minn. 
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- it is reasonable that a closer correlation could not readily be obtained 

fomavery limited mmber oftests. 
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3 right purchase ofr radar equipment that has been designed specifically for — 


R. BOWMAN,! M. ASCE. 
_ comprehensive progress report on the status f the use of radar as a means 
for making quantitative precipitation forecasts. He also should be thanked , 
his sparing use of the involved terminology that has become the language o 
_ the radar meteorologist. Although radar meteorology is still an infant a 7 
‘years, a tremendous amount of research has been undertaken in this | field; call a 
lengthy bibliography appended to this paper merely highlights some ‘of the 
more significant developments in the subject field of endeavor. Owing to the ss 


a nature of its original presentation, »| the paper is necessarily b brief, and and many ye 


has come e to recognize r many of the catenitel uses, as well as the limitations, 


of radar as a meteorological instrument. The greatest deterrent to the out- 


meteorological use probably is its first cost. By way of illustration, the cost 
of each high-power, 10-centimeter, weather radar unit currently being 
in service by the United States Weather Bureau amounts to about 
$100,000, exclusive « of the cost of associated structures In n 
of such high costs, most of the weather radar units currently used 
institutions other than Federal governmental agencies, are modified 
_ Government- -surplus units that were designed originally for military duty, but 
- which can be purchased for a fraction of their original cost. In military appli- aii 
a we precipitation echoes are not desirable because they tend to mask or _ a 
"obscure the echoes received from the primary targets; for this reason, the 
‘ circuitry in military radar units is designed to enhance the echoes from the Ls 


_ desired targets and to suppress precipitation echoes. Modification of a mili- _ 
q tary ‘radar unit is largely confined to the alteration of the video presentation be 4 


receiver or transmitter circuits could involve the redesign of the entire ——- 
unit ‘Thus the advantage of low first cost of surplus radar r equipment (the — 


= 


Sc in the receiver section, inasmuch as any significant change in basic — 


- author’ s “discarded equipment”) £ frequently is offset by compromises in per- 

. a author’ s complaint concerning “equipment of improper wavelength” is _ 

justifiable, , particularly where it concerns 3-cm radar applied to quantitative 
work, It has been determined experimentally that for showers of 16 = 
tensity, (18) the signal | attenuation | at 7 cm is tiiantail 60 times greater than 


Proc. 1901, January, 


| 

— 


10 cm. . At greater ‘intensities this is even more = 
- Numerous studies of field situations of various types lend support to the grow- 


ae The major deterrent to the more extensive use of 10-cm units again is the 
‘first cost. ‘Dimensions of the larger components of a radar unit are about | 
ie: proportional to the wavelength, so that the volume occupied by a complete 10- 
cm unit may be several times seme than the space occupied by a 3- -cm unit: t 
From the foregoing discussion, it might appear that 3- -centimeter radar 
has little value as a meteorological instrument; this is not the case. Although 
the 3-cm band may not be the optimum band for quantitative prediction of ees 
precipitation at long range, it has considerable utility in determining the areal 
extent and vertical development of potentially severe storms. This has been | 
_ demonstrated conclusively by the author and his associates in Illinois, and 
through its adoption as a key unit in the tornado tracking networks of the _ acl 
Southwest. Some of the 3- “cm units originally designed for military airborne | 
use, ‘such as the AN/APS- 4, are so compact that they } are readily converted to to 
small mobile units and easily moved into relatively inaccessible areas; ‘such — 
7 unit has been put to effective | use by Decker(19) and — in the re of ah 2 


parisons drawn in rainfall depths are compared with 
4 contours observed over a short period during each hour. _ This would imply _ F 
_ that the rainfall intensities, as well as the distribution of heavy oor orga 

centers, as indicated by radar echoes, ‘were assumed to remain constant dhe’ 


ie . throughout each hourly period. On the contrary, rainfall is notoriously varia- | 


ble in practically all types of storms. The writer suggests that an integrated 
series of time-lapse photographs, similar to those associated with the studies 
of mesosystems undertaken by the University of Chicago,(20) might represent 
more closely the ever-changing storm patterns. 
es The development of the pulse > integrator (illustrated by Fig. 4) represents — 
a a giant step toward accurate correlation of echo and rainfall intensities. i Pa 
st principal disadvantage, however, is that the time required for gating and inte- 
limits its effectiveness to relatively small areas. 
_ The nature of the received power and rainfall intensity ; traces s of Fig. 5 
wa considerable promise for accurate correlation. The errors inherent in 
comparisons of instantaneous values appear to be of a compensating nature, - 
which suggests the possible improvement of accuracy by a time-integration — 
method. _ Values of signal intensity could be converted to represent db values" 
- above a dbm cutoff value. In an approximate analysis of the left-hand group 
of curves, the writer obtained very | close correlation by shifting the time- 
integrated rainfall curve 2 minutes to the left, and by integrating the signal — z 


to the significance ¢ of the cutoff value, but he the that tit 
“may ay be related to signal attenuation. | 


HY 
— 
20. 
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4 wa 
paper illustrates several of the methods by which correlation of radar 
echoes with observed precipitation has been attempted; all illustrate the need 
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Rockney, v.D., “The WSR- 57 Radar,’ Proc. Seventh Weather Radar ‘Seventh Weather 


Contf., Amer. ‘Met. , Soc. and Univ. Miami (Fia.), 1908. 


‘Purcell, E. M., “The Radar Equation, Radar System Engineering (L. N. 
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19. "Decker, F. W. and E Pike M. 4 Radar Investigations of Orographic 
Precipitation, ” Pate. Seventh Weather Radar Conf. 
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by E. E. Pariset and B. 


E M. ASCE and B. Rhone’s very interesting 
"paper has brought up a number of the use of low head 


“a _ Among the largest radial gates built ; is the 148 ft wide by 30 ft high gate 


built for Donzere Mondragon (France), 
_ We had the opportunity to test, on scale models, several spillways using — 


: ie gates, so we are happy to add a few “stones” to the —* built t by 

About discharge characteristics, we generally use the fourth method: __ 

Qs C D Ly 2g H, as any formula is difficult to use for very large gate — 


ings when the discharge depends more of the crest profile than of the gat —_ 
opening, and as for low openings, true orifice flow occurs. = B49 
a The value of C depends not only of 6 (angle formed by the tangent to | = 2 
gate lip and the tangent to the crest curve at the nearest point of the crest cere 
curve), but also of the crest shape which governs to a certain extent the 
4 pressure on the spillway at the gate outlet. _ The thinner the crest profile sa 
comparison with the “datum shape ” the more likely some depression is to oc- 
cur and the higher will be the coefficient of discharge, as shown on graph * 


(Brazil). At full’ gate opening the coefficient of discharge, , using: gall 
‘This profile was tested and used for the surelevation of the dam. ol ies 
- a profile much thinner than the “datum shape”, graph No. 1, was used not — 
only to obtain a higher discharge and to save on the volume of concrete, but © 
5 mainly t to allow the new profile of the surelevated dam to be tangent with the 
: old one, , in order to avoid the water nappe from jumping free of the down- = 
stream face of the dam and falling out of the — ‘bucket at the foot of the 
must be noted that varies also with the value H, a greater head in- 


‘creasing the depression and the coefficient of discharge. 


_ All our tests have confirmed that if the gate is located just on the on 


is sufficient to the pressure for ali gate openings higher or equal 
- to the m minimum pressure occurring at full gate openings. yeas 
Proc. Paper 1935, February, 1959, by Thomas J. Rhone. 
1. Director, LaSalle Hydr. Lab., Montreal, Canada. 


Research , LaSalle Hydr. Lab., Montreal, Canada. 
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DISCHARGE COEFFICIENT Head to center ot gote 


ing downstream the gate sealof Edgardde SouzaDam. = ~~ 
e ‘The problem is exactly the same if if a vertical lift gate | is used instead of a 
#3, radial one. _ We had the opportunity to test such an an arrangement for the Hart- 
aaa To determine the maximum allowable pressure on the profile ofa ene, 
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where: p/ p/w 


water ‘vapor pressure 
velocity head at the point of minimum 


could be transformed locally into additional negative pr behind 


constructed very smoothly so that no irregularity will be 


_— if the spillway will operate under conditions giving these low | 
= 1. 0 no danger of cavitation 


, Cavitation on weir crests, PP. Cc meeting. 
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LA IONS HIP Ss BETWEEN PIPE RESISTANCE FORMULAS? 


by Peter Ackers 


PETER ACKERS. Rs. —The a author 8 ’s historical review ew of pipe resistance 


at present. Unfortunately its complexity has prevented this 
coming into general use, there = two main obstacles | to its adoption in 


I (a) It cannot be expressed asa a grid of ‘straight lines on a simple design 


_ Ife, v, R- von s are quoted, V may be calculated directly, but in the 


design case of known and available atr ial 

The Moody3) and Rouse(4) non-dimensional plots of the equation, , though 
"clarifying the picture, were only a partial solution of the problem facing the ! 
engineer. Neither Lamont’ (5) addition of nomograms nor Bonnington’ s(6) 

_ transposition of the variables succeeded in reducing the application to the : rs 

7 simple form required by the hydraulic engineer, i.e., , design charts oo 
in terms of conduit size, hydraulic gradient, velocity discharge. 
However, studies of the propiem at the Hydraulics Research Station have pro-- 


_ The physical non-dimensional parameters, Nr, f and e/ D are replaced by 
non- dimensional parameters which varia- 


V, D and S, as follows: 
_ : fi. a 


Ss, in terms of these new 


— 
4 
in the United Kingdom, which drew the similar conclusion that the _ 
= 
q 
at 
: 
te 
mij 
— 
Principal Scientific Officer, Hydrs. Research Station, Wallingford, 


dimensional, it retains some physical significance and applies | 
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any system of units. = stric form, occur explicitly, so that 
The engineering variables, in parametric form, oc 
— for a given class of pipe Suid, 
— (iii) A resistance diagram applicable to any or row ameters, Fig. 1. 
— on the basis of these dimensionless parameter ed 
nt 
|. 
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7 


For any roughness and kinematic viscosity, this non- Fa; 

_ dimensional chart can readily be converted to a 1 design chart in terms Cs 

z of V, R and S, merely by tracing the e appropriate acctien ant re- labelling - 
_ (v) The addition of a grid of discharge lines provides a | complete direct so- i 


lution to flow problems in which any two of the four basic variables, 


> 


There now seems no excuse for replacing the accepted resistance equation 


_ by exponential approximations, except in certain limited fields of hydraulic x iD 
_ computation, such as the pipe-network example cited by the author. _ However, © 
‘in ref. 7, the writer demonstrated how this can be done by partial anageilae 


central values of R “ine D) and S. 


the between a and may be: written 


where R, is the roughness ni mabe, wil It follows Eq. (8) that 
the oft-quoted criterion dimensional 
is true e zero, i.e., it is applicable only to smooth- con- 
4 ditions. This also results from Eq. (9) by putting R, equal to zero, whilst ae 


The assumption implicit in the author’s method of deducing exponential 
"approximations from the Moody diagram is that relative roughness is constant. 


, f _ However, for any given class of commercial pipe, it is ; the absolute roughness . 


_ which should be regarded as constant, and this is an important distinction if _ 
an equation covering various pipe diameters is required. To this extent, the 
‘writer method to general than the author’ 

REFERENCES a. 

ee a “Resistance « of fluids flowing in channels and pipes”, Hydraulics 
Research ch Paper No. 1, London, .S.0., 1958. 


Cc. F. Colebrook, “Turbulent in pipes, with particular to the 
transition region between smooth and rough laws”, Instn civ. 


& 
— 
| 
a 
= 
ay 
i 
_ 


OF OS On 08 @« S$ @ 


sec for 


i sad jaa} 


Vide 


WN 


0-002 
Oe 


having coment or ciphelt cocting 


tree 
men 
joan 


usecs) for PIPES flowing full 


construction rough form, rovgh 


que 

ing 


oot % 


GE 


7 
“ 
273 


DISCHARGE Q 


— AS AN 
: 
— 
— 


“Friction f for pipe | Sow", A. Mech. 
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